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Abstract 



Quantum Dots: Coulomb Blockade, 
Mesoscopic Fluctuations, and Qubit Decoherence 

The continuous minituarization of integrated circuits is going to affect the un- 
derlying physics of the future computers. This new physics first came into play as 
the effect of Coulomb blockade in electron transport through small conducting is- 
lands. Then, as the size of the island L continued to shrink further, the quantum 
phase coherence length became larger than L leading to mesoscopic fluctuations - 
fluctuations of the island's quantum mechanical properties upon small external per- 
turbations. Quantum coherence of the mesoscopic systems is essential for building 
reliable quantum computer. Unfortunately, one can not completely isolate the system 
from the environment and its coupling to the environment inevitably leads to the loss 
of coherence or decoherence. All these effects are to be thoroughly investigated as 
the potential of the future applications is enormous. 

In this thesis I find an analytic expression for the conductance of a single electron 
transistor in the regime when temperature, level spacing, and charging energy of an 
island are all of the same order. I also study the correction to the spacing between 
Coulomb blockade peaks due to finite dot-lead tunnel couplings. 1 find analytic 
expressions for both correction to the spacing averaged over mesoscopic fluctuations 
and the rms of the correction fluctuations. 

In the second part of the thesis I discuss the feasibility of quantum dot based 
spin- and charge-qubits. Firstly, I study the effect of mesoscopic fluctuations on the 
magnitude of errors that can occur in exchange operations on quantum dot spin- 
qubits. Mid-size double quantum dots, with an odd number of electrons in the range 
of a few tens in each dot, are investigated through the constant interaction model 
using realistic parameters. It is found that the number of independent parameters 
per dot that one should tune depends on the configuration and ranges from one to 
four. Then, I study decoherence of a quantum dot charge qubit due to coupling to 
piezoelectric acoustic phonons in the Born-Markov approximation. After including 
appropriate form factors, I find that phonon decoherence rates are one to two orders of 
magnitude weaker than was previously predicted. My results suggest that mechanisms 
other than phonon decoherence play a more significant role in current experimental 
setups. 



m 



Acknowledgements 

First and foremost, I would like to thank my adviser Prof. Harold U. Baranger 
for his thoughtful guidance on every aspect of my research and for the funding of my 
work all these years. I am also grateful to Prof. Konstantin A. Matveev for investing 
so much time in my education and his guidance on my first research project at Duke. 
Prof. Eduardo R. Mucciolo has been a wonderful collaborator and friend. He gave 
me the confidence that I can publish in Physical Review B. Without him this thesis 
would have been much thinner. 

I am grateful to Profs. Berndt Miiller and M. Ronen Plesser for helping me to 
secure TA positions when I needed them most. 

I am grateful to the members of my Ph.D. committee: Profs. Shailesh Chan- 
drasekharan, Albert M. Chang, Gleb Finkelstein, and Weitao Yang for finding time 
to be on my committee and valuable comments on the manuscript. 

During my years at Duke, I benefited from stimulating discussions with Alexey 
Bezryadin, Alexander M. Finkelstein, Martina Hentschel, Alexei Kaminski, Eduardo 
Novais, Stephen W. Teitsworth, Denis Ullmo, Gonzalo Usaj, and Frank K. Wilhelm. 

I would like to thank my friends: Sven Rinke, Alex Makarovski, Kostya Sabourov, 
Anand Priyadarshee, Sung Ha Park, Ji-Woo Lee, Martina Hentschel, Gonzalo Usaj, 
Eduardo Novais, Ribhu Kaul, and Oleg Tretiakov. Thank you, guys, for being there 
for me even when I did not ask. It is always fun to be around you! 

Finally, a special thank you goes to my parents, grandparents, and my wife Elena 
whose patience helped me to get the job done. And, certainly, life makes much more 
sense because my daughter Tanya is around. 



IV 



Contents 



Abstract iii 

Acknowledgments iv 

Coulomb Blockade in Quantum Dots 1 

1 Introduction to Quantum Dot Physics 1 

1.1 Overview 1 

1.2 2D Lateral Quantum Dots 1 

1.3 Constant Interaction Model and Single-Particle Hamiltonian 3 

1.4 Constant Exchange And Interaction Model 4 

1.5 3D Quantum Dots 6 

2 Coulomb Blockade Oscillations of Conductance at Finite Energy 
Level Spacing in a Quantum Dot 8 

2.1 Overview 8 

2.2 Introduction 8 

2.3 The Model 10 

2.4 Linear Conductance in the Spinless Case 12 

2.5 Application to Tunneling Through Quantum Hall Edge States in a 
Quantum Dot 18 

2.6 Linear Conductance in the Spin-| Case 21 

2.7 Conclusions 32 

3 Coulomb Blockade Peak Spacings: Interplay of Spin and Dot-Lead 
Coupling 35 



3.1 Overview 35 

3.2 Introduction 36 

3.3 The Hamiltonian 40 

3.4 Plan of the Calculation 41 

3.5 Doublet-Triplet-Doublet Spin Sequence: Calculation of the Spacing 
Between Peaks 42 

3.5.1 Zeroth Order: Isolated Quantum Dot 43 

3.5.2 Second Order: Contribution From Virtual Processes 43 

3.6 Ensemble Averaged Correction to the Peak Spacing 50 

3.7 RMS of The Correction to Peak Spacing due to Mesoscopic Fluctuations 54 

3.8 Conclusions 58 

Quantum Computation with Quantum Dots: Mesoscopic Fluctuations 

and Decoherence 59 

4 Introduction to Quantum Computation with Quantum Dots 60 

4.1 Overview 60 

4.2 Introduction to Quantum Computation 60 

4.3 Mathematical Basis of Quantum Computations 61 

4.3.1 Quantum Bits and Entangled States 61 

4.3.2 Steps in Quantum Computation and Quantum Gates 62 

4.3.3 Quantum Parallelism 63 

4.3.4 Important Remarks 64 

4.4 DiVincenzo's Criteria of Physical Implementation 65 

4.5 Quantum Computing Proposals 65 

4.5.1 Liquid State NMR 65 

4.5.2 Spintronic Quantum Dot Proposal by Loss and DiVincenzo . . 66 

5 Spin Qubits in Multi-Electron Quantum Dots 69 

5.1 Overview 69 

5.2 Introduction 69 

5.3 Model System 71 

5.4 Errors in Exchange Operations 73 

vi 



5.5 Results 75 

5.5.1 Mesoscopic Effects 76 

5.5.2 Relevance for Real Quantum Dots 79 

5.6 Conclusions 80 

6 Time Evolution of the Reduced Density Matrix 81 

6.1 Introduction 81 

6.2 Formulation of the Problem 82 

6.3 General Solution 82 

6.4 r^o = and Born Approximations 84 

6.5 Markov Approximation on Top: Redfield Equation 85 

7 Phonon Decoherence of a Double Quantum Dot Charge Qubit 87 

7.1 Overview 87 

7.2 Introduction 87 

7.3 Model System 89 

7.4 Decay of Charge Oscillations 93 

7.5 Bias Pulsing 97 

7.6 Conclusions 100 

8 Conclusion 101 

8.1 General Overview 101 

8.2 Main Conclusions of My Work 101 

8.3 Possible Directions of Future Research 104 

Appendices 105 

A Occupation Numbers in the Canonical Ensemble 106 

B On the Electron-Phonon Coupling in GaAs quantum dots 108 

B.l Introduction 108 

B.2 Deformation Potential Contribution to the 

Electron-Phonon Coupling 108 

vii 



B.3 Piezoelectric Contribution to the 

Electron-Phonon Coupling 110 

B.4 Combining Two Mechanisms Ill 

C Derivation of Eqs. (7.18)-(7.20) 112 

Bibliography 114 

Biography 124 



vm 



To my grandmother Vorozhtsova, Tatyana Ivanovna (1917-1987) 



IX 



Chapter 1 

Introduction to Quantum Dot 
Physics 

1.1 Overview 

Quantum dot research [1, 2, 3, 4] has developed into an exciting branch of meso- 
scopic physics. Many novel phenomena were observed in transport measurements 
through quantum dots: Coulomb blockade [5], even-odd asymmetry in Coulomb 
blockade peak spacings [6, 7], and Kondo effect [8, 9] to name just a few. The field 
has kept researchers busy for about twenty years now and still continues to surprise 
us. 

By no means can 1 provide a detailed introduction to the entire field in this 
thesis (even restricting myself to quantum dot physics alone) and I do not think it is 
necessary as there are quite a few nice review papers available [5, 10, 11, 12]. Instead, 
I give just enough introductory material on 2D lateral and 3D quantum dots so that 
the reader can jump into the chapters where the original results of my research are 
presented. 

1.2 2D Lateral Quantum Dots 

Recent advances in materials science made possible the fabrication of small con- 
ducting devices known as quantum dots. In particular, in 2D lateral quantum dots 
from one to several thousand electrons are confined to a spatial region whose linear 
size is from about 40 nm to l/xm [2, 13]. These quantum dots are typically made 
by (i) forming a two-dimensional electron gas on the interface of semiconductor het- 





GaAs AL Ga^.j^As 



^1|J^n 



Figure 1.1: 2D lateral quantum dot. Left panel: Schematic picture of the quan- 
tum dot setup. Two-dimensional electron gas is formed on the the AlGaAs/GaAs 
heterostructure interface. Metal electrodes (or gates) are deposited on top of the 
heterostructure. Negative voltages are applied to the electrodes to repel the electron 
gas underneath - bright regions contain no electrons. Then, bias voltage is applied 
between source and drain leads (see wires coming out of the sample) and the current 
through the quantum dot is measured. Right panel: Micrograph of a real quantum 
dot [4]. Bright regions correspond to the deposited metal electrodes. Arrows show 
the flow of electrons. Copyright by C M. Marcus, from Ref. [4]. 



erostructure and (ii) applying electrostatic potential to the metal surface electrodes to 
further confine the electrons to a small region (quantum dot) in the interface plane [4], 
see Fig. 1.1. 

The transport properties of a quantum dot can be measured by coupling it to 
leads and passing current through the dot. The electron's phase is preserved over 
distances that are large compared with the size of the system (quantum coherence), 
giving rise to new phenomena not observed in macroscopic conductors. 

The coupling between a quantum dot and its leads can be experimentally con- 
trolled. In an open dot, the coupling is strong and the movement of electrons across 
dot-lead junctions is classically allowed. However, when the point contacts are pinched 
off, effective barriers are formed and conduction occurs only by tunneling. In these 
almost-isolated or closed quantum dots, the charge of the dot is quantized and the 
dot's low-lying energy levels are discrete with their widths smaller than the spacing 
between them, see Chapter 2. 

The advantage of these artificial systems is that their transport properties are 
readily measured and all the parameters - the strength of the dot-lead tunnel cou- 



plings, the number of electrons in the dot, and the dot's size and shape - are under 
experimental control. 

To observe quantization of the quantum dot charge, two conditions have to be 
satisfied. Firstly, the barriers must be high enough so that the transmission is small. 
This gives the following condition for the conductance: G <^ e^//i, that is, the dot 
must be almost or completely isolated. Secondly, the temperature must be low enough 
so that the effects of charge quantization are not washed out. The quantum dot's 
ability to hold charge is classically described by its average capacitance C. Since 
the energy required to add one electron is approximately e^/C, we find the following 
condition: T <^ e^/C. 

The tunneling of an electron onto the dot is normally blocked by the classical 
Coulomb repulsion with the electrons already in the dot; hence, the conductance 
is very small. This phenomenon is known as the Coulomb blockade. However, by 
changing the voltage of the back-gate Vg one can compensate for this repulsion and, 
at the appropriate value of Vg, the charge of the dot can fluctuate between n and n + 1 
electrons leading to a maximum in the conductance. Thus, one can observe Coulomb- 
blockade oscillations of the conductance as a function of the back-gate voltage, see 
Fig. 2.4. At sufficiently low temperatures these oscillations turn into sharp peaks that 
are spaced almost uniformly in Vg. Their separation is approximately equal to the 
charging energy e^/C. 

1.3 Constant Interaction Model and Single- Particle 
Hamiltonian 

Electron-electron interactions in a quantum dot in the Coulomb-blockade regime 
are conventionally described by the constant interaction (CI) model [5]. In this model 
the Hamiltonian of the system is given by a sum of two terms: (i) the electrostatic 
charging energy, which depends only on the total number of electrons in the dot n 
and (ii) the Hamiltonian of free quasiparticles: 

Hci = Ecn"^ + '^£c,cl^Caa, (1-1) 

aa 

where Ec = e^/2C and Caa is the annihilation operator of a quasiparticle (electron) 
on orbital level a with spin a. 

Single-particle dynamics inside real 2D lateral quantum dots with more than 40 
electrons has no particular symmetry due to irregular boundaries (chaotic quantum 
dot) or the presence of impurities (disordered quantum dot). In both of these cases, 
free quasiparticles inside the dot can be described by random matrix theory (RMT): 
The Hamiltonian is chosen "at random" except for its fundamental space-time sym- 
metry [14, 15, 16, 17, 18, 10]. Random matrix theory is applicable if the dimensionless 



conductance of the dot g is large: g = Et/Se ^ 1, where Et is the Thouless energy 
and Se is the mean level spacing in the quantum dot. For a ballistic quantum dot 
Et ~ hvp/L, where vp is the Fermi velocity of the electrons and L is the linear size of 
the dot. A large value of g indicates that the dot can be treated as a good conductor. 

The spacings between conductance peaks contain two contributions as well. The 
first one, due to the charging energy, does not fiuctuate much. The second contri- 
bution is proportional to the spacing between discrete energy levels in the quantum 
dot. This term does fluctuate and obeys the Wigner-Dyson statistics. Spin degener- 
acy of each one-particle energy level in the quantum dot leads to the even-odd parity 
effect: the second contribution appears only when we promote an electron to the next 
orbital. 



1.4 Constant Exchange And Interaction Model 

More careful treatment shows that the interactions between electrons in a quan- 
tum dot should be correctly described by the "universal" Hamiltonian [19, 12]. 

In the basis of eigenfunctions {(pa} of the free-electron Hamiltonian, 



-J-v^ + f/W 



,(r) = £«0„(r), (1.2) 

the two-particle interaction takes the form 

where U{r) is the random potential determined by the shape of the quantum dot and 
o"i and o"2 are the spin indices for the fermionic operators. The generic matrix element 
of the interaction is 

Hap-yS = / rfri(ir2y(ri - r2)0a(ri)0;3(r2)0*(r2)0Kri)- (1-4) 

The matrix elements of the interaction Hamiltonian have a hierarchical structure. 
Only a few of these elements are large and universal, whereas the majority of them 
are proportional to the inverse dimensionless conductance 1/g {1/g = 27r//cp'L for 
a 2D quantum dot) and, therefore, small [19, 12]. As a result, the Hamiltonian 
[Eq. (1.3)] can be broken in two pieces: 

Hint = -f^ini + -f^i„/ • (1-5) 

The first term here is universal - it does not depend on the quantum dot geometry and 
does not fiuctuate from sample to sample (for samples differing only by realization 



of disorder). The second term in Eq. (1.5) does fluctuate but it is of order 5e/ g and, 
hence, small. This term only weakly affect the low-energy {E < Et) properties of the 
system. 

The form of the universal term in Eq. (1.5) can be established using the require- 
ment of compatibility of this term with the RMT [19, 12]. Since the random matrix 
distribution is invariant with respect to an arbitrary rotation of the basis, the opera- 
tor /j}J may include only the operators which are invariant under such rotations. In 
the absence of the spin-orbit interaction, there are three such operators: 

^ = ^^ciaCaa (1-6) 

- the total number of electrons, 

Q.Ul(72 

- the total electron spin of the dot, and the operator 

CatCai, (1-8) 



E 



which corresponds to the interaction in the Cooper channel. 

Gauge invariance requires that only the product T'^T of the operators T and 
T'^ may enter the Hamiltonian. At the same time SU(2) symmetry dictates that 
the Hamiltonian may depend only on S^ and not on the separate spin components. 
Taking into account that the initial interaction Hamiltonian [Eq. (1.4)] is proportional 
to c^ we find the "universal" Hamiltonian: 

Hf;l = Ecn' - J.S^ + J,f tf , (1.9) 

where Ec is the redefined value of the charging energy [19, 7] and Js > is the ex- 
change interaction constant. The constants in this Hamiltonian are model-dependent. 
The first two terms represent the dependence of the energy on the total number of 
electrons and the total spin, respectively. Because both the total charge and the total 
spin commute with the free-electron Hamiltonian, these two terms do not have any 
dynamics for a closed dot. The situation changes as one couples the dot to the leads. 

The third term vanishes in the Gaussian Unitary (GUE) random matrix ensemble. 
GUE corresponds to the absence of time-reversal invariance, or placing the quantum 
dot into an external magnetic field. One can say that the Cooper channel is suppressed 
by a weak magnetic field (it is sufficient to thread a unit quantum flux $o through 
the cross-section of the dot). 



Thus, in the absence of superconducting correlations, the universal part of the 
interaction Hamiltonian consists of two parts: 

VcEi = Ecn" - J,S\ (1.10) 

This is the so-called constant exchange and interaction (CEI) model [20, 7]. Its dom- 
inant part depends on the QD charge number n - the corresponding energy scale 
Ec = e^/2C is related to the capacitance of the QD, C, and exceeds parametri- 
cally the mean level spacing 6e. The second part depends on the total spin S - the 
corresponding energy scale Jg is less than 6e. 

If the level spacings did not fluctuate, then the smallness of Jg would automat- 
ically imply that the spin of the QD can only take the values of (if n is even) or 
I (if n is odd). Fluctuations in the level spacings may lead to a violation of this 
periodicity [21, 22, 19]. However, the Stoner criterion; Jg < 6e, guarantees that the 
total QD spin is not macroscopically large; that is, it does not scale with the QD 
volume. 

In the presence of the back-gate electrode capacitively coupled to the QD, the 
CEI model and free quasiparticle Hamiltonian become 

ij(°) = Y^^kUka + Ec{h- Nf - JsS^ (1.11) 

ka 

where M = CgVg/e is the dimensionless back-gate voltage and Cg is the dot-backgate 
capacitance. 



1.5 3D Quantum Dots 

A 3D quantum dot is just a metal nanograin. There are numerous methods to 
synthesize them [23] . At the time of writing, metal nanoparticles of down to about 
1 nm in diameter are commercially available. 

To form a single electron transistor one has to trap a nanograin between two 
leads. To create a narrow gap between two conductors, the electromigration technique 
is often implemented [24]. In this technique one takes a conductor of a small cross- 
section and gradually increases the voltage and, hence, the current through it until 
cracking occurs. Using this method, one can create very narrow nanometer size gaps. 
To trap a nanoparticle the potential difference between two leads is applied. Then the 
nanoparticle gets polarized and attracted to the region of high electric field. This is 
called the electrostatic trapping technique. Thus, the nanoparticle is attached to two 
leads via oxide tunnel barriers. As an example, in the recent experiment by Bolotin 
and coworkers an ultra-small gold nanograin, 5 nm in diameter, was incorporated into 
a gap between two leads. Thus, a single-electron transistor was formed, and Coulomb 
blockade oscillations were observed for more than ten charge states of the grain. 



Although these metal nanograins are similar to semiconductor quantum dots, 
there are a number of important differences [11]. (i) Metals have much higher den- 
sities of states than semiconductors; hence, they require much smaller sample sizes 
(less than 10 nm) before discrete energy levels in the QD become resolvable. The 
ratio Ec/5e is usually larger for nanograins; therefore, mesoscopic fluctuations have 
significantly less impact on their quantum properties, (ii) For nanograins, the tun- 
nel barriers to the leads are formed by insulating oxide layers. Therefore, they are 
insensitive to applied voltages, whereas for 2D quantum dots they can be tuned by 
changing voltages on the electrodes. 



Chapter 2 

Coulomb Blockade Oscillations of 
Conductance at Finite Energy 
Level Spacing in a Quantum Dot 

2.1 Overview 

In this chapter we find an analytical expression for the conductance of a single 
electron transistor in the regime when temperature, level spacing, and charging energy 
of a grain are all of the same order. We consider the model of equidistant energy levels 
in a grain in the sequential tunneling approximation. In the case of spinless electrons 
our theory describes transport through a dot in the quantum Hall regime. In the case 
of spin-i electrons we analyze the line shape of a peak, shift in the position of the 
peak's maximum as a function of temperature, and the values of the conductance in 
the odd and even valleys. 

2.2 Introduction 

Recent progress in mesoscopic fabrication techniques has made possible not only 
the creation of more sophisticated devices but also greater control over their proper- 
ties. Electron systems confined to small space regions, quantum dots, and especially 
their transport properties have been studied extensively for the last decade [11, 5]. In 
particular, an individual ultra-small metallic grain of radius less than 5 nm was at- 
tached to two leads via oxide tunnel barriers, thus forming a single electron transistor 
(SET) [11]. Applying bias voltage, V, between two leads allows one to study transport 




Figure 2.1: Scheme of the Coulomb blockade setup. 



properties of the system, Fig. 2.1. Alternatively, a SET can be formed by depleting 
two-dimensional electron gas at the interface of GaAs/AlGaAs heterostructure by 
applying negative voltages to the metallic surface gates [5]. 

In this chapter we will assume that the bias voltage is infinitesimally small, V -^ 
0. This corresponds to the linear response regime. In order to tunnel onto the 
quantum dot, an electron in the left lead has to overcome a charging energy, Ec = 
e^/2C, where e > is the elementary charge; C is the capacitance of the quantum 
dot. If T <^ Ec then conductance through the system is exponentially suppressed. 
This phenomenon is called the Coulomb blockade. However if we apply a voltage, 
Vg, to the additional gate capacitively coupled to the dot, the Coulomb blockade can 
be lifted. Indeed, changing Vg one can shift the position of energy minimum so that 
energies of the quantum dot with iVg and A'^g + 1 electrons will become equal and 
an electron can freely jump from the left lead onto the dot and then jump out into 
the other lead. Thus, current event has occurred and a peak in the conductance, G, 
corresponding to this gate voltage is observed. By changing the gate voltage one can 
observe an oscillation of the conductance or Coulomb blockade oscillations. 

One-particle energy levels in the quantum dot, {Ei}, are given by the solution of 
the Schrodinger equation in the quantum dot's potential. The mean spacing between 
these energy levels is 6E. The conventional assumption that Ec ^ SE is not valid 
in the case of sufficiently small dots. In fact, in the recent experiments [25, 26, 27], 
where a Cqq molecule has acted as a quantum dot, the level spacing is of order 
charging energy. Experiment [25, 26] was performed at T = 4.2K as well as at room 
temperature. In other experiments [28, 29] with quantum dot formed by depleting 
2DEG [28] and ropes of carbon nanotubes acting as a quantum dot [29], charging 
energy is only three times larger than the spacing 6E. 
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Though Coulomb blockade oscillations have been studied in a number of impor- 
tant limiting cases [30, 31, 32, 33], the problem in the case when values of Ec-, SE, 
and T are all of the same order has not been theoretically addressed. Let us note 
that energy levels of the quantum dots are random and obey Wigner-Dyson statistics 
with the fluctuation of order of their mean [18]. Nonetheless to go as far as possible 
in the analytical treatment of the problem we have to assume that energy levels in 
the quantum dot are equidistant. In this chapter we derive an analytical expression 
for the linear conductance, G = I/V\v-,o, in the case of spinless as well as spin-| 
fermions. 

In Section 2.3 we describe our model and the assumptions involved. We write the 
model assuming spin-i fermions. In Section 2.4 we consider the linear conductance in 
the case of spinless fermions. We obtain an analytical expression for the conductance 
and analyze its limiting cases. In Section 2.5 we consider one possible application of 
the Section 2.4 results, namely tunneling through the edge states in a quantum dot 
placed into a strong magnetic fleld. In Section 2.6 the linear conductance as well as its 
properties in the case of spin-| fermions is considered. In Section 2.7 we summarize 
our flndings. 



2.3 The Model 

Hamiltonian of the system in question is 

H = Hi + Ha + f. (2.1) 

Here, the flrst term is the Hamiltonian of noninteracting electrons in the left and 
right leads: 

Hi = Y, Ekclcka + J2 EpcU^'pa, (2.2) 

ka per 

where a continuum of states in each lead, \ka), \pa) is assumed; E^, Ep and Cka, 
Cpa are the energies and electron annihilation operators in the left and right leads, 
respectively; a stands for the z-component of spin. The chemical potentials of the 
leads, jj, ^ Ec,SE,T, are shifted according to the bias voltage, V, applied. Fig. 2.2. 
We will assume that leads are in thermal equilibrium at temperature T and, thus, 
occupied according to the Fermi-Dirac distribution. 

The second term in Eq. (2.1) is the Hamiltonian of the quantum dot: 

Hd = J2E^clc,^ + U, (2.3) 



tcr 



where flrst term is the kinetic energy of electrons in the quantum dot: {Ei} is a 
discrete set of the quantum dot's energy levels; Qo-'s are the annihilation operators. 
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^i 



Figure 2.2: Electrostatic potential energy along a line through the tunnel junctions. 



The second term, tJ describes the electron-electron interaction in the quantum dot. 
We adopt the simplest model for the interaction, namely, the constant interaction 
model. In this model the Coulomb interaction of the electrons depends only on the 
total number of electrons in the quantum dot: 



U{N) = EcN' - eVeN, (2.4) 

where N = ^^^ Ci^^i^ ~ Ni is the total number of excess electrons; Ni is the total 
number of positively charged ions. The second term is the contribution from external 
charges. They are supplied by the ionized donors and the gate: Ve = Vd + aVg, where 
a is a function of the capacitance matrix elements of the system. Thus, Ve can be 
varied continuously by changing gate voltage, Vg. U{N) can be rewritten as 



U{N) = Ec{N - Ngf + Const, 

where Ng = eVe/2Ec is the dimensionless gate voltage. 

The third term in Eq. (2.1) is the tunneling Hamiltonian: 



(2.5) 



X^ [tkicl^Ci^ + h.c. 



/ J y^pi^^pcr^icr ' il'C'j 



(2.6) 



kiu 



ptcr 



where tki and tpi are matrix elements of tunneling into the left and right leads, re- 
spectively. 

We assume that the dot is weakly coupled to the leads; that is, the conduc- 
tances of the dot-lead junctions are small: G'''^ <^ ^'^/h, where h is Planck's constant. 
Equivalently, the widths of the quantum dot's energy levels contributing to the con- 
ductance, Fj = r' -|- r[, must be small compared to spacing between them: Fj <^ 6E. 
This, together with Fj ^ T assumption, allows us to characterize the state of the dot 
by a set of occupation numbers, {njo-} [33]. 
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2.4 Linear Conductance in the Spinless Case 

The model formulated above has been studied by Beenakker in the sequential 
tunneling approximation [33]; that is, conservation of energy was assumed in each 
tunneling process, and cotunneling was neglected. Therefore, to find the stationary 
current, kinetic equation considerations can be applied. In the linear response regime 
an analytical formula for the conductance has been obtained. In the case of spinless 
fermions [33]: 

x[l-nF{Ei-fi + U{N)-U{N-l)], (2.7) 



where 



and 



r^ = 27r J] \tk^\^6 [E, -Ek + U{N) - U{N - 1)] 






r[ = 2n}^\tpi\'5 [Ei -Ep + UiN) - U{N - 1)] 

V 

are widths of the quantum dot's level i associated with tunneling into the left and 
right leads, respectively; Peq{Ne) is the equilibrium probability that the quantum 
dot contains iVg electrons; Feq{Ei\Ne) is the occupation number of level i given that 
the dot contains N^ electrons; npiE) is the Fermi-Dirac distribution; and fi is the 
chemical potential in the leads. 

The quantity Fgg{Ei\Ng) in (2.7) is the most non-trivial one to calculate. It is the 
occupation number of the level i in the canonical ensemble (A'^e is fixed). In the limit 
SE/T -^ 0, Feq{Ei\Ne) bccomcs a Fermi-Dirac distribution with the appropriately 
chosen chemical potential: fl = {Eq + Ei)/2, where Eq corresponds to the energy of 
the last occupied energy level at T = 0, Fig. 2.3(a); Ei corresponds to the energy of 
the first empty energy level at T = 0. In the opposite limit 6E/T -^ cxd, the Fermi- 
Dirac distribution with /i = {Eq + Ei)/2 apparently breaks down: the occupation 
number of level j = 1, for example, see Fig. 2.3a, is rii = e^^^^'^, not e^*^^/^-^ as the 
Fermi-Dirac distribution would predict [33]. 

The occupation number in question is 



FeqiE,\N,) = ^^-L_J2Pe,i{ni})6^^,i6^^,^ 



{ni} 

mA^.)^e-^S^-5„.,5^.,j,„,, (2.^ 

{ni] 



where Peg({^/}) is the equilibrium probability of the |{^/}) state of the quantum dot; 
(3 = 1/T; and the detailed definition of F{Ne) will follow. The reason for writing this 
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j = 3 
j = 2 
j = l 
j = 



j = l 
j = 
j = -l 
J = -2 



a) 



b) 



Figure 2.3: a) occupation numbers in tlie canonical ensemble for the equidistant 
energy levels at T = 0; b) mapping sum over i onto sum over j {N = 2 case is 
shown) . 



equation is to show that analytical calculation of the occupation numbers is hardly 
possible for arbitrary quantum dot's energy level structure, {Ei}. 

The only way to overcome this difficulty is to assume that energy levels in the 
quantum dot are equidistant. Then one can use the bosonization technique [34] (see 
Appendix A) to find the exact analytical expression for the occupation numbers in 
the canonical ensemble. It was done by Denton, Muhlschlegel, and Scalapino [35]: 



Feq{Ei\N, 



n~ 



m=l 



.nni-lg-5[m^- 



-(2i 



'l)m\^ 



where 



Ei 



J 



C 1 
- — ^ + - 
5E 2 



integer. 



(2.9) 



(2.10) 



6E/T = 6 = const, ( is the energy corresponding to the highest occupied energy level 
at T = plus SE/2, Fig. 2.3(a). This quantity ( is somewhat similar to the chemical 
potential of a dot, though, strictly speaking, the chemical potential is not well-defined 
for a dot in the canonical ensemble. The difference, C — /x, is a linear function of the 
gate voltage. Therefore, by properly adjusting "zero" value of the gate voltage it can 
be put to zero. Hereinafter, we assume that ( — fi = 0. 

To calculate the conductance we also need to know Peg{Ne), the probability that 
the dot, in thermodynamic equilibrium with the reservoirs, contains N^ electrons. It 
can be calculated in the grand canonical ensemble: 



Peg{N, 



7 ' 



(2.11) 
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where Z^ = Yllv =0^"'^'^^^^^ ^^ ^^^ grand partition function; '^{N^) is the thermody- 
namic potential of the quantum dot. It can be expressed via free energy of the dot's 
internal degrees of freedom, F{Ne): 



hence, 



where 



V9(iV,) = F(iVe) - fiN, + U{N), (2.12) 



P,,{N,) = ^e-^^(^)Z(iVe), (2.13) 



is the partition function in the canonical ensemble. In the last expression the sum is 
taken over all possible states, |{nj}) of the quantum dot. To calculate Z{Ne) exphcitly 
we need to assume that energy levels of the dot are equidistant: 

Z(iV,) = e-(^-^>)'^/2^e.c, (2.14) 

where iVg = Ni corresponds to the equilibrium number of the excess electrons {( = fi); 
Zg^c is the partition function of the thermal excitations. Now, let us substitute 
Eq. (2.14) in Eq. (2.13): 

= _]_^-miN)+N^SE/2]^ (^2.15) 

where 

00 
D',= J2 e-^l^^^^+^'-^^/^i^ (2.16) 

Af=-oo 

Here, we have extended one limit of the sum to infinity since the Fermi energy is the 
largest energy scale of the problem. 

Thus, we are well-equipped to calculate the conductance in the case of equidistant 
energy levels in the quantum dot at arbitrary ratio S. The widths of energy levels, 
r' and r[, in the quantum dot are energy dependent, random quantities. Let us 
assume that quantum dot is weakly coupled to the leads via multichannel tunnel 
junctions: G'''' = CfN^h <^ e^/^j where N^h is the number of channels; Gf is 
the conductance of one channel. Experimentally, this situation corresponds to the 
metallic grain coupled to the leads via oxide tunnel barriers [11]. This setup allows 
one to decrease fluctuations of the energy levels' widths, T\ and r[, by a factor of 
yJNch- We also assume that the widths are slowly changing functions of the energy, Ei. 
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Then, T] and T^ can be considered constants and evaluated at the chemical potential: 
r' ^ r' = r'; r^ ^ rj^ = r*". There is a simple relation between these widths and 
conductances of the corresponding junctions. In the case of spinless fermions: 

Ti = ^SE, r = ^5E. (2.17) 

Let us substitute Eq. (2.15) in Eq. (2.7): 

G = ^|^^E-"'^'^"^""^'"'^'E^-(^^i^e) 

N i 

x{l-nF[Ei-fi + U{N) - U{N - 1)]} . (2.18) 

To take advantage of the expression for occupation numbers, Eq. (2.9), we need to 
map the sum over i onto the sum over j. As illustrated in Fig. 2.3, the mapping rule 
depends on the total number of excess electrons, N (compare with Eq. (2.10) written 
for N = 0): 

E,-fi = E,-C' + {C-C) = (j-l)sE + N6E, 

where (' = C{^) is the energy of the highest occupied energy level in the dot with 
N excess electrons at T = plus SE/2. We will also use the following identities: 

5E . / 5E\ . 



U{N) + —N' = [Ec + ^)N- eVeN 
SE\ / , , 1^ ^ 



Ec + — 1 ( iV - Ao + - 1 + Ci, (2.19) 

where Ci = - {eV^f /2{2Ec + 5E)- 

eVe I , , 

has been chosen so that Aq = corresponds to the maximum of the conductance 
peak; and 

U{N) - U{N - 1) = (2A^ - l)Ec - eVe 
= (2iV - l)Ec - {2Ec + 5E)(Ao - 1/2). (2.21) 

Substituting these results in Eq. (2.18), we obtain 

Goo ~ D,^ Z.e-^-^-2(Jv-Ao)so + i' ^^■^^) 
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Figure 2.4: Coulomb blockade oscillations of the conductance as a function of the 
dimensionless gate voltage, Aq at 5E = Eq- Curves are plotted for different temper- 
atures: T/Ec = T/6E = 1/6 = 0.2, 0.35, 0.5, 0.7, 1, 2, 5. 



where 



Do = exp{pCi)D'^ = J2 



g-eo(JV-Ao+|; 



(2.23) 



N 



eo = PiEc + 5E/2); G^ = G^G''/{G^ + G') is the classical, Ec, 6E <^T <^ fi, limit 
of the conductance. We have also used the identity: 1 — np{E) = (1 + e"^^)""*^. 
Eq. (2.22) is the general expression for the linear conductance in the spinless case for 
equidistant energy levels in the quantum dot at arbitrary values of Ec, SE and T. 

One can immediately prove the following properties of the conductance, Eq. (2.22). 
First of all, G(Ao) = ^(Ao + M), where M is an integer. In the gate voltage units, 
AVe = {2Ec + SE)/e is a period of the conductance oscillations. This property re- 
flects symmetry with respect to adding (removing) an electron to the quantum dot. 
Secondly, due to the electron- hole symmetry, conductance is an even function of Aq: 
G{Ao) = Gi-Ao). 

The linear conductance, Eq. (2.22), as a function of the dimensionless gate voltage 
Aq at 6E = Eq is plotted in Fig. 2.4 for different temperatures. At low tempera- 
tures there are sharp Coulomb blockade peaks. At high temperatures, T ^ Eq, SE, 
Coulomb blockade is lifted and small oscillations of the conductance can be observed. 
These oscillations are slightly non-sinusoidal and given by the following asymptotic 
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Figure 2.5: The exact line sliape of tlie conductance peak, Eq. (2.26), plotted at 
1/5 = T/5E = 0.4 (solid curve). Dotted curve corresponds to blindly applying 
classical regime formula (2.27) at T/6E = 0.4. Dashed curve corresponds to blindly 
applying formula (2.28) at T/SE = 0.4. The argument of the plot is linear function 
of the gate voltage, 0o- 



formula: 



GjAo) - G(Ao) ^ 2^3/^6 



^eo/4 



Gr 



-7r2/eo 



cos(27rAo 



+ e-2"'/"" cos(47rAo) + ©(e"^"'/"") 



(2.24) 



where G(Ao) is the average value of the conductance. The second term in Eq. (2.24) 
is due to inherently non-sinusoidal nature of the conductance oscillations, see Fig. 2.4. 
To derive this expression one can use Poisson's summation formula. 

To study the line shape of a separate peak let us consider the limit of large 
charging energy: Ec ^ T, 6E or, equivalently, Eq ^ 1,6 in the dimensionless units. 
In Eq. (2.22) only A^ = —1,0 terms in the sum over A^ give substantial contribution to 
the conductance near Aq = 0; all other terms are exponentially suppressed. Besides, 
the sum over j at the A^ = — 1 is 0(e~^^") and, therefore, can also be neglected. 
Hence, line shape of the conductance peak at Aq = is given by 






6 



-2£oA 



A^ I 



+ e- 



-i5+2eoAo * 



(2.25) 
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It is more instructive to rewrite this equation as follows: 

(2.26) 



^(00) ^ S v^ nj 

Goo ~ 1 + e-'^^'^o ^ 1 + e-J''^+^^'^o ' 
i 

^(0). T.(0) 



where (po = Ve — Ve ; Ve is chosen so that 0o = corresponds to center of the 
conductance peak. In the classical regime [30], T ^ 6E, line shape of the conductance 
peak is given by 

G'(0o) ^ /^e0o 
Goo 2sinh(/3e0o)' ^ ' ' 

In the opposite limit of 6E ^ T: 

^(00) 5 



Goo 2[l + cosh(/3e(/ 



(2.28) 



The exact line shape of the conductance peak, Eq. (2.26), at T/6E = 0.4 is shown in 
Fig. 2.5. On the same figure we also plotted two conductance peaks in the limiting 
cases, Eqs. (2.27) and (2.28), out of their validity region at T/6E = 0.4. Nevertheless, 
it is interesting that the exact conductance peak is higher than both of the limiting 
cases peaks. The peak's height is given by 

G(0) _ 5 ^ n, _ / 1/2, T » 6E ^^.29) 



Goo 2 ^^ 1 + e-i^ 1 5/4:, 5E > T > T 



3 



Temperature dependence of the conductance peak's height was numerically calculated 
in the Ref. [33], see Fig. 2 there. 



2.5 Application to Tunneling Through Quantum 
Hall Edge States in a Quantum Dot 

Formulas for the linear conductance in the case of equidistant energy levels in 
a dot and spinless fermions derived in Section 2.4 can be applied to a number of 
physical problems. 

Let us consider, for example, a quantum dot formed by confining a two-dimensional 
electron gas by a circularly symmetric electrostatic potential, U{r). We assume that 
U{r) is zero at the origin and takes large value at r = i?, where R is the radius of the 
dot. Fig. 2.6(a). Let us apply a strong magnetic field, B, perpendicular to the plane 
of the dot. This situation corresponds to the quantum Hall regime and was reviewed 
in Ref. [36]. 

To solve the one-electron Schrodinger equation in this geometry it is convenient 
to choose the symmetric gauge. Then, angular momentum is, clearly, an integral 
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a) 



b) 



Figure 2.6: Geometry and spectrum of the quantum dot states: a) symmetric gauge 
eigenstates, |m), including edge state are shown schematically; b) energy spectrum of 
the eigenstates with different angular momenta, two lowest Landau levels are shown. 



of motion. In each Landau level, ra, states with larger angular momentum, \m), 
are localized further from the origin, near a circle with radius Rm = /_f/\/2(m + 1), 
where Ih = \/hc/eB is the magnetic length, c is the speed of light. The presence 
of the confinement potential leads to an increase in energy for the symmetric gauge 
eigenstates with Rm of order or larger than R (or, m of order or larger than mp, 
see Fig. 2.6(b)). For the states with Rm ~ i? an electron is influenced by both the 
electric field of the boundary, E{R) = U'{R)/e, and strong, perpendicular to the 
electric, magnetic field. Thus, near the edge electron executes rapid cyclotron orbits 
centered on a point that slowly drifts in the direction of E x B, that is, along the 
boundary. Thus, Quantum Hall edge states are formed. Fig. 2.6(a). It is important 
to notice that in this closed geometry electron system has only one edge. In this 
consideration we also assume that Ih "C R. 

For simplicity let us consider the case when only the zeroth Landau level crosses 
the chemical potential, that is, there is only one type of edge states. This corresponds 
to a sufficiently strong magnetic field so that filling factor, u is equal to 1, Fig. 2.6(b). 

Now, we are ready to consider transport through this type of quantum dot in the 
strong magnetic field. Let us weakly couple it to two leads and apply an infinitesimally 
small bias voltage between them. An electron from the left lead can now tunnel into 
dot's edge state and then tunnel into the right lead as illustrated by dashed lines in 
Fig. 2.6(a). 

The energy spectrum of edge states can be linearized as follows 






dS 

om 



[ra — rap) 



(2.30) 



mp 



20 



Thus, energy levels of the edge states are equally spaced with the spacing [36] 



8E 



dSr, 



dm 



d£„ 



mp 



dRr, 



dR„ 



R 



dm 



nip 



eEiWf 



(2.31) 



This fact makes formulas derived in Section 2.4 applicable to this problem. Essential 
assumption here is that dispersion curve, Fig. 2.6(b) is almost linear in the range of 
angular momentums: \m, — m^pl ^ Am,, where Am, = m.a.x{l , T / 6 E) . 

In the case at hand, spacing, 6E is inversely proportional to the size of a dot just 
like charging energy, Ec- Hence, their ratio does not depend on the size of a dot and 
is given by 



6E 
Ec 



tE{R) 



a 



B 



(2.32) 



where e is the dielectric constant of the media around the interface; a = e^ /he is the 
fine structure constant. Therefore, in this case oscillations of the conductance given 
by Eq. (2.22) are determined by only one parameter 5 = 6E/T. 

In conclusion, let us consider the case of an arbitrary shaped quantum dot. In 
this case, m is just the index of an edge state and no longer associated with the 
angular momentum. The phase along the boundary for the m,-th edge state is 



Or. 



(A/tJj '^TfT, \ "-' I 1 



(2.33) 



where km{x) is the corresponding wave vector, x parametrizes the boundary, and L 
is its length. The phase difference between two consecutive edge states is 



9. 



m-\-l ^m 



9m = 2n 



dx[km+l{x) - km{x)] , 



(2.34) 



where km+i{x) — km{x) = {Sm+i — £m) /hv{x), v{x) = cE{x)/B is a drift speed along 
the boundary. Then, the spacing between edge states' energy levels is 



SE = 27ch 



dx 

v{x) 



27re/| 



dx 



(2.35) 



Though the electric field E[x) at the boundary slightly changes as one goes from one 
edge state to the other, this effect is small and we neglect it. Therefore, the energy 
levels of the edge states are equidistant with the spacing given by Eq. (2.35). 

In the case of the circularly symmetric quantum dot, E{x) is constant, and one 
can easily perform the integration in Eq. (2.35). This leads to the previously obtained 
expression for the level spacing, Eq. (2.31). 
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2.6 Linear Conductance in the Spin-^ Case 

Formula (2.7) for the linear conductance in the spinless case can be easily gener- 
alized to the spin-i case by counting each energy level twice [33]: 

X [1 - npiEi -fi + U{N) - U{N - 1))], (2.36) 

where Feq{Ei-i\Ne) is the occupation number of the quantum dot's energy level i with 
a spin-up electron, (i, t) in the canonical ensemble: number of electrons in the dot, 
A'^e, is fixed. 

As in the spinless case, to carry out analytical consideration we have to assume 
that energy levels in the quantum dot are equally spaced. Presence of the spin 
degeneracy makes the calculations more complicated. 

First of all, let us find the occupation number Feq{Ei-^\Ne). Let us consider spin- 
up and spin-down electron subsystems. Ground state energy of the system is 

E, = {N^ + Nf) ^ = [(AT, + N,f + (AT, - N,f] ^, 

where N„ = {Ne)„ — Ni/2 is the number of excess electrons in the spin-a subsystem; 
A^j is chosen even. N^ + Nj^ = N is the total number of excess electrons in the 
quantum dot; (A^| — N^)/2 = S^ is 2;-component of the total electron spin. Using 
these identities, one can find that 

Eg = jN^6E + Sl6E. (2.37) 

While Sz is subjected to the thermodynamic fiuctuations, A^ is fixed. 

The occupation number in question, Feq{Ei^\Ne) = n^i, is known if, in addition 
to A^, the z-component of the total spin, 5*2, is fixed. In the case of an even number 
of electrons, parameter ^' for the spin-up electron subsystem is equal to C + SzSE, 
given Sz, see Figs. 2.7(a) and 2.7(b). Occupation numbers in the spin-up subsystem 
at fixed Sz are given by Eq. (2.9) with the appropriately chosen parameter (': fij^Sz- 

However, 2;-component of the total spin, Sz, is not fixed but subjected to the 
thermodynamic fiuctuations. Therefore, to find the occupation numbers, n^p we 
have to account for all possible values of Sz'- 



Ne/2 

ev 

5z=-Afe/2 



nT,= Y. P{Sz)n,.s., (2.38) 
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Figure 2.7: Parameter (' of the spin- up electron subsystem for even number of 
electrons: a) at Sz = 0, b) at arbitrary integer Sz {Sz = 1 case is shown); for odd 
number of electrons: c) at Sz = 1/2, d) at arbitrary half-integer Sz {Sz = 3/2 case is 
shown) . 
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where 

-Sl5 

nS.) = -^, (2.39) 

S,=~Ne/2 

is the probabihty that ^-component of the total spin of the quantum dot is equal to 
Sz- Substituting Eqs. (2.39) and (2.9) into Eq. (2.38) we obtain 

Ne/2 






ev 



/2 

m=i Y2 e-'5'i'5 

Since number of electrons in the quantum dot is even, Sz may take only integer values. 
Therefore, 



%T 



oo 
ev 



y^ tl\m-l^-[^+U-i)m]5 S,=~oo 



m=i 2~] e 

Sz=~00 



-S?5 



where we extended limits of the sum over Sz to infinities since the Fermi energy is 
the largest energy scale in the problem. Separating m = 2r — 1 and m = 2r parts of 
the sum, where r is a positive integer, we obtain final expression for the occupation 
numbers in the case of even number of electrons: 



oo oo 

3\ 



nT, = A{6) J2 e-(^-5)('-+2.--|)5 _ ^ ^-r(r+2j-l)S^ (2.40) 



ev 

r=l r=l 



where 



(s-h)h 



E e- 

AiS) = '-^^ . (2.41) 



-oo 



In two limiting cases 

r 2e-^/4 [1 + (e--^)] , 5E^T 

A{5) = \ 

I 1 - 4e-^'/^ + O ('e-^-VA ^ T ^ 6E 

Analytical expression for the high-temperature limit of A[5) can be obtained using 
Poisson's summation formula. 
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One can easily prove the following properties of the occupation numbers n^t valid 
at arbitrary temperature: 



rf;\ = l-n-^-p (2.42) 

e^' Vf + e"-'' V4.. = A{6)e^/\ (2.43) 



They are valid due to the electron-hole symmetry and similar to the following proper- 
ties of the Fermi-Dirac distribution: npiE) = 1— nF{—E) and e^^nF{E)+e~^^np{—E) 
= 1. 

Similarly, one can find occupation numbers in the case of odd number of electrons 
in the quantum dot, Ng. Energy level which contains one electron at T = will be 
referred to as j = level. In this case electron-hole symmetry corresponds to j — > — j 
transformation. Parameter (' of the spin-up electron subsystem at a given Sz is equal 
to C + {Sz - l/2)6E, see Figs. 2.7(c) and 2.7(d). Therefore, 







< = 


Ne/2 

■ >: 


P{Sz) 

12 


^■- 


-(^-1; 


00 
m=l 


-1)™- 


1 ^(^ 


■+jm\S S, 


N^/2 

E e 


-(5. 


= -f)'^ 


N,/2 

E 

Sz = -Ne/ 


e" 
2 


-s'is 



(2.44) 



Since number of electrons in the quantum dot is odd, Sz may take only half-integer 
values. Separating odd and even parts of the sum over m, we obtain: 



'^ " A{5) 



-.00 00 

i^^e-(^-^)K2^-5)'5_^e-^(^+2,)5_ (2.45) 



r=l r=l 



Property of the electron- hole symmetry reads as follows: n°i = 1 — n"^,^ 

It turns out that there exists simple relation between rf^ and n°'^ occupation 
numbers: 

nf^ = A(5)e~(^"3) V^-^. (2.46) 

This property is the analog of upiE) = e~^^np{—E) one of the Fermi-Dirac distri- 
bution. It will allow us to get rid of n°'^ occupation numbers in the final expression 
for the conductance. 

Now we are in a position to find the probability that a dot, in thermodynamic 
equilibrium with the reservoirs, contains A'^e electrons, Peg(^e)- Eq- (2.13) written 
for the spinless case is still applicable if we keep in mind that energy levels in the 
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quantum dot are doubly degenerate. Partition function of the dot's internal degrees 
of freedom in the canonical ensemble is 

= E E Z{Nl)Z{N^)^^^^^U^,, (2.47) 

nJ=on}=o 

where 

Z{Nl) = J2^''^^^'~'^'"^NlEn. (2-48) 

{rii} 

is the partition function of the spin-up electron subsystem in the canonical ensemble. 
Mathematically, expression for Z{Nl) is identical to the one for Z{Ne) in the spinless 
case. Thus, one can directly apply the result obtained previously, Eq. (2.14): 

Z{Nl) = e-(^'-^'/2)25/2^^^^^ (2.49) 

where Nj = Ni/2 is the equilibrium number of electrons in the spin-up subsystem. 
Similar result is valid for the partition function of the spin-down electron subsystem, 
Z{N^). Substituting these results in Eq. (2.47) we obtain: 



Ne N, 



nJ=on^=o 
The exponent can be simplified as follows 



+N^,N^- 



hence, 



■7fAT\ ^^ ^<= 



N^ Ne/2 

y^ g-(ivI-jv./2)25 ^ ^^N^5/A y^ 

iVj=0 S=-NJ2 



^-NHIA y^ g-(iVl-JV,/2)25 ^ g-iV2V4 V^ ^-sH 



where in the second equality we took advantage of the delta symbol. Sum in the last 
line is taken over integer values of s if N^^ is even or half-integer values of s if A^e is 
odd. Limits of the sum over s can be extended to infinities since we assume that 
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fi ^ T. According to Eq. (2.13) probability that quantum dot, in thermodynamic 
equihbrium with the reservoirs, contains N^ electrons is 






5 



^^ ^^ s=~N./2 

g-/3f/(Af) J- g-s^5 
s=-Ne/2 



N=ev s N=od s 



where U{N) = U{N) + N'^6E/4; and we used the fact that A^ and N^ have the same 
parity since Aj is chosen even. Therefore, sums over N = ev and N = od are taken 
over A^ = 0, ±2, ±4, . . . and A^ = ±1, ±3, . . . values, respectively. At this point in 
the calculation we need to specify whether the total number of electrons in the dot 
is even or odd: 

n,l^Vej-| (^jj,yl^(^s^^-mN)^ A^ is odd ' ^ ' 

where 

D'= J2 e'"''^''^ + A{6) Y, e-'"'^^'^^. (2.51) 

N=even N=odd 

Now we are prepared to calculate the conductance, Eq. (2.36), in the case of 
the equidistant double degenerate energy levels in the dot at an arbitrary 5E/T and 
Ec/5E ratios. Similarly to the consideration in the spinless case we assume that 
quantum dot is weakly coupled to the leads via multichannel tunnel junctions, and 
tunneling widths of the energy levels in the quantum dot, T\ and r[, are slowly 
changing functions of the energy, Ei. Then, these tunneling widths can be considered 
constants and evaluated at the chemical potential: V\ ^ V^^ = F'; F^ ^ FJ^ = F**. 
Furthermore, they can be expressed via conductances of the corresponding junctions: 



, _ hG' 6E _ hG^ SE 

^ ~ ~72 V ^ ~ ~72 V \^.^^) 



There is an additional factor of 1/2 here compared to the spinless case, Eq. (2.17), 
due to the double degeneracy of each energy level in the quantum dot. First of all, 
let us break the sum over A'e in Eq. (2.36) in two parts: A^ = even and A^ = orfrf, and 
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f I j-Q j = 2 j = l 

f I T'T'' J-Q f I j--i 
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a) b) c) 

Figure 2.8: Mapping of the sum over i onto the sum over j for a) A^ = 0; b) even 

A^ (N = —4 case is shown); c) odd N {N = —3 case is shown). 



apply Eqs. (2.52) and (2.50): 

Gr 



G" + G"- D' 

oo 

^ e-'^^(^) ^Fe,(E,tl^e) {1 -np[E,-fi + U{N) - U{N - 1)]} 



X 

\N=even i=l 



+A{5) J2 e-'^^(^) ^Fe,(E,t|iVe) {1 -nF[E,-f, + U{N) - U{N - 1)]} .(2.53) 



N=odd i=l 



To take advantage of the occupation numbers we derived, Eqs. (2.40) and (2.45), we 
need to map each of the sums over i onto the sum over j. The first sum over i in 
Eq. (2.53) is taken at even number of excess electrons, see Figs. 2.8(a) and 2.8(b), 
hence 



F,,{Ei^\N,) = n 



ev 



E,-i, = E,-C+iC-C) = (j-l)5E + ^NSE. 

Remember that by properly choosing "zero" of the gate voltage we put ( = fi. The 
second sum over i is taken at odd number of excess electrons, see Figs. 2.8(a) and 
2.8fc), therefore 



Feg{Ei^\Ne) = n' 



od 



E,-fi = E,-C+{C-0 = j6E+^N6E. 
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We will also use the following identities: 

U{N)= fEc+^-^)N'-eVeN= ^E^ + ^"j (n-A + ^] + C2, 
where C2 = -(eV;)V(4^c + SE); 

is the dimensionless gate voltage, A = corresponds to a position of the conductance 
peak; and 

U{N) - U{N - 1) = (2iV - l)Ec - eVe = (2iV - l)Ec - 2{Ec + ^)(A - ^). 
Substituting these results in Eq. (2.53) we obtain 

G{A) ^ ^ Sr^ ^-e{N-A+l/2f 

n Z^ 

n 

1 1 



\2 



Goo D 

N=even 



e 

'=even 

CO 



X 5^ n: 



ev 
j = -oo 



,(2.55) 



-(j-l/4)<5-2(Af-A)e _^ ]_ g-(j-l/4)5+2(Af+l-A)e _^ I 
j = _oo 

where 

D = exp(/5C2)D'= J] e-^(^-^+i/2)^ + A((5) 5^ e-^(^-^+i/2)'; (2.56) 



N=even N=odd 



e^p(Ec + ^-pj; (2.57) 

Goo = G"G'7(G' + G"') is the high-temperature, Ec,SE <^ T <^ fi, limit of the 
conductance. To eliminate n°'^ from the final expression we used useful property of 
the occupation numbers given by Eq. (2.46). 

Formula (2.55) is the main result of this chapter. It is the analytical expression for 
the linear conductance in the spin-| case for equidistant energy levels in the quantum 
dot. One can use Eq. (2.55) to plot Coulomb blockade oscillations of the conductance 
as a function of the dimensionless gate voltage. A, at arbitrary values of Ec, SE 
and T. Particularly, when all energy scales are of the same order: Eq ~ SE ~ T, 
numerical calculation is a breeze. Fig. 2.9. 

One can immediately notice the following properties of the linear conductance. 
First of all, G(A) = G(A + 2M), where M is an integer. In other words, {4:Ec + SE)/e 
is the conductance period in the gate voltage units. This property reflects symmetry 
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Figure 2.9: Linear conductance oscillations as a function of the dimensionless gate 
voltage, A, at bE = Eq- Curves are plotted for different temperatures: T / Ec = 
T/6E = 1/6 = 0.15, 0.3, 0.5, 0.7, 1, 2, 5 using Eq. (2.55). 



with respect to adding (removing) two electrons to (from) the quantum dot. Secondly, 
conductance is a symmetric function with respect to the center of a valley. A' = 
M+1/2, where M is an integer. That is, G(A- A') = G(A'- A). This is a reflection 
of the electron-hole symmetry. These properties of the conductance oscillations are 
not generic. They are valid due to the assumption of equally spaced energy levels in 
a quantum dot. 

At high temperatures, T ^ Ec, SE, conductance peaks overlap and their max- 
imums become almost equidistant. Fig. 2.9. As a result, instead of separate peaks, 
the conductance in this limit has oscillatory behavior. 

Let us find temperature dependence of the conductance in the valleys. In the 
sequential tunneling approximation the conductance in the valleys decays exponen- 
tially as T ^ 0, Fig. 2.10. At low temperature number of electrons in a dot in the 
valleys is almost quantized. We will call the valley "odd" ( "even" ) if it corresponds 
to odd (even) number of electrons in the dot. We find that at any temperature the 
conductance in the odd valley is larger than that in the even one. Fig. 2.10. This 
feature is robust with respect to the distribution of energy levels in a quantum dot. 

However, it is important to mention that at low temperatures, T < Tj„, where 
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Figure 2.10: Temperature dependence of the conductance in tlie odd, G'(— l/2)/G'oo 
and even, G(l/2)/G'oo valleys at bE = Ec- Upper curve corresponds to the conduc- 
tance in the odd valley. 



cotunneling [9] will dominate sequential tunneling contribution to the conductance 
in the valleys. Therefore, temperature dependence of the conductance in the valleys. 
Fig. 2.10, is valid only for the temperatures T > Tj„. 

Let us analyze the limit of large charging energy: Eq ^ T, SE or, equivalently, 
e ^ 1,6 in the dimensionless units. In this limit, two adjacent peaks in the conduc- 
tance have exponentially small, ~ e~^'^^'^ , overlap with each other. Thus, it makes 
perfect sense to study the line shape of a separate peak. Let us determine hne shape 
of the conductance peak near A = 0. In the numerator of Eq. (2.55) only the N = 
term in the sum over A^ survives; moreover, at A^ = second term in square brackets 
is 0(e~^^). In the denominator, only the A^ = —1,0 terms matter. Hence, the line 
shape of the conductance peak near A = at arbitrary 6E/T ratio is given by 



g(0) 



6 



A{5)e 



-I3ecf> 



,^r 



"-it 



-(j-l/4)5+/3e<^' 



(2.59) 



where we used the following identity: 

2£A = f3e (K - V}°^) = f3e(f). 

Clearly, = corresponds to A = 0. In the classical regime, T ^ 6E, the line shape 
of the conductance peak is given by [30] 

G{(j)) jSecj) 



Goo 2 sinh(/9e(; 



(2.60) 
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Figure 2.11: The exact line sliape of tlie conductance peak, Eq. (2.59), plotted 
at 1/(5 = T/5E = 0.4 (solid curve). Dotted curve corresponds to blindly applying 
classical regime formula, Eq. (2.60), at T/SE = 0.4. Dashed curve corresponds to 
blindly applying formula (2.61) at T/6E = 0.4. The argument of the plot is linear 
function of the gate voltage, 0. 



Formally, this equation is identical to that of the spinless case, Eq. (2.27). Nonethe- 
less, the values of Goo are different in these two cases by a factor of 2. This is due 
to spin degeneracy of each energy level in the spin-i case, compare Eqs. (2.17) and 
(2.52). In the limit of 5^ > T [31]: 



G(0) 



Goo 3 + 2y/2cosh{pe(f)+\S-^\n2)' 
This peak has its maximum at 

e<t>LT = e0„ (T < 6E) = -^6E + -^T, 



(2.61) 



(2.62) 



and is symmetric with respect to this value: G(0 — 0lt) = G{<Plt — 4>)- The exact 
line shape of the conductance peak, Eq. (2.59), at T/6E = 0.4 and two limiting 
cases conductance peaks, Eqs. (2.60) and (2.61), plotted out of their vahdity region 
at T/6E = 0.4 are shown in Fig. 2.11. As in the spinless case, the exact conductance 
peak is higher than both of the limiting cases peaks. 

The position of the peak's maximum, 0m = 4>m{T), is shifted to the left from its 
high temperature limit, (pcL = (pm(T ^ 6E) = 0. It is determined by the equation for 
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Figure 2.12: Temperature dependence of the dimensionless peak's maximum posi- 
tion, e(f)m{T)/SE. In the low temperature hmit: e(f)m{0)/SE = —1/4 according to 
Eq. (2.62). 



(pm'- G'i'Pm) = 0, where G(0) is given by Eq. (2.59). The dimensionless position of the 
peak's maximum, e(j)rn{T)/5E, as a function of the temperature, T /5E, is numerically 
plotted in Fig. 2.12. 

The conductance peak height is Gmax = G{<Pm)- In the limiting cases: 



Gr 



G{<Pr, 



(jTr 



Czn 



1/2, T>5E 
{3-2^2)6, 6E:^T:^Ti 



(2.63) 



Peak's height as a function of the temperature, T/6E, can be plotted numerically. 
Fig. 2.13. 



2.7 Conclusions 



We have studied Coulomb blockade oscillations of the linear conductance through 
a quantum dot weakly coupled to the leads via multichannel tunnel junctions in the 
sequential tunneling approximation. To obtain analytical results we have assumed 
that the energy levels in the dot are equally spaced. The electron-electron interaction 
in a quantum dot has been described by the constant interaction model; though, 
thermal excitations with all possible spins have been taken into account. 

The linear conductance in the spinless case is given by Eq. (2.22). It is valid 
at arbitrary values of Ec, SE and T. The line shape of an individual conductance 
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Figure 2.13: Height of tlie conductance peak, Gmax/Goo-, as a function of tlie tem- 
perature, T jbE. Dotted and daslied curves correspond to two limiting cases peak 
heights, Eq. (2.63). 



peak at arbitrary ratio b = 6E/T is given by Eq. (2.26). Exact conductance peak is 
higher than both of the hmiting cases peaks at any gate vohage as is illustrated in 
Fig 2.5. An analytical expression for the height of the conductance peak at any ratio 
S is obtained, Eq. (2.29). 

In Section 2.5 we applied the spinless case theory result to the problem of the 
transport via a dot in the quantum Hall regime. Energy levels in a dot in this case 
are equidistant with the spacing given by Eq. (2.35). 

Linear conductance in the case of spin-| electrons at arbitrary values of Eq, SE, 
and T is given by Eq. (2.55). In particular, this equation allows one to plot the 
conductance oscillations in the regime when the charging energy, level spacing in the 
dot, and the temperature are all of the same order. Fig. 2.9. We find that the period 
of Coulomb blockade oscillations is doubled compared to the model with a continuous 
electronic spectrum in the dot. Equation (2.55) is the main result of the chapter. 

We also find that conductance in the odd valley is larger than that in the even 
one at any temperature. Fig. 2.10. The difference between conductances has the 
largest value at T ^ 0.3 6E (at 6E = Ec). The sign of the difference is the same 
as for the quantum dot in the Kondo regime [9]. Kondo effect takes place at very 
low temperatures, T < Tk <^ Tj„, where Tk is the Kondo temperature, and leads 
to the logarithmic enhancement of the conductance in the odd valleys [12]. Our 
consideration shows that even-odd asymmetry exists at much higher temperatures. 
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Line shape of the conductance peak is given by Eq. (2.59). As in the spinless 
case, the conductance peak is higher than both of the hmiting cases peaks at any 
gate voltage, Fig 2.11. As we increase the temperature peaks' maximums shift and 
become more equidistant. Fig 2.12. The peak's height as a function of the temperature 
is calculated numerically and plotted in Fig. 2.13. 

Though we have found physical system which has equidistant energy levels in 
the spinless case, see Section 2.5, we are not aware of any such system in the spin-i 
case. In the case of a chaotic quantum dot Wigner-Dyson model gives a fairly good 
approximation for the distribution of the energy levels of the dot. If we had assumed 
Wigner-Dyson distribution of the quantum dot's energy levels then we would have 
had to give up the hope of finding a solution. It goes back to the very difficult 
problem of finding occupation numbers of the dot's energy levels in the canonical 
ensemble. The only way to solve it is to assume that energy levels in the quantum 
dot are equally spaced. Then one can use the bosonization technique to find the 
occupation numbers. Assumption of the equidistant energy levels is in line with the 
level repulsion property of the Wigner-Dyson distribution. Therefore, the analytical 
consideration of this reasonably simplified model, in our opinion, is a significant step 
forward in the solution of the general problem. 

Though we do not expect our quantitative results to precisely describe a quantum 
dot with random energy levels, they certainly give correct order of magnitude for the 
conductance oscillations and their generic features. 



Chapter 3 

Coulomb Blockade Peak Spacings: 
Interplay of Spin and Dot-Lead 
Coupling 

3.1 Overview 



For Coulomb blockade peaks in the linear conductance of a quantum dot, we 
study the correction to the spacing between the peaks due to dot-lead coupling. 
This coupling can affect measurements in which Coulomb blockade phenomena are 
used as a tool to probe the energy level structure of quantum dots. The electron- 
electron interactions in the quantum dot are described by the constant exchange and 
interaction (CEI) model while the single-particle properties are described by random 
matrix theory. We find analytic expressions for both the average and rms mesoscopic 
fluctuation of the correction. For a realistic value of the exchange interaction constant 
Js, the ensemble average correction to the peak spacing is two to three times smaller 
than that at Js = 0. As a function of Js, the average correction to the peak spacing 
for an even valley decreases monotonically, nonetheless staying positive. The rms 
fluctuation is of the same order as the average and weakly depends on Jg. For a 
small fraction of quantum dots in the ensemble, therefore, the correction to the peak 
spacing for the even valley is negative. The correction to the spacing in the odd 
valleys is opposite in sign to that in the even valleys and equal in magnitude. These 
results are robust with respect to choice of the random matrix ensemble or change in 
parameters such as charging energy, mean level spacing, or temperature. 

The work in this chapter was done in collaboration with Harold U. Baranger. 
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3.2 Introduction 



Progress in nanoscale fabrication techniques has made possible not only the cre- 
ation of more sophisticated devices but also greater control over their properties. 
Electron systems confined to small regions - quantum dots (QD) - and especially 
their transport properties have been studied extensively for the last decade [5, 11]. 
One of the most popular devices is a lateral quantum dot, formed by depleting the 
two-dimensional electron gas (2DEG) at the interface of a semiconductor heterostruc- 
ture. By appropriately tuning negative potentials on the metal surface gates, one can 
control the QD size, the number of electrons n it contains, as well as the tunnel barrier 
heights between the QD and the large 2DEG regions, which act as leads. Applying 
bias voltage V between these leads allows one to study transport properties of a single 
electron transistor (SET), Fig. 3.1(a) [5]. 

We study properties of the conductance G through a QD in the linear response 
regime. We assume that the dot is weakly coupled to the leads: GL,/?*^e^//i, where 
Gl,r are the conductances of the dot-lead tunnel barriers, e > is the elementary 
charge, and h is Planck's constant. 

To tunnel onto the quantum dot, an electron in the left lead has to overcome a 
charging energy Eq = e^/2C, where C is the capacitance of the QD, a phenomenon 
called the Coulomb blockade. However, if we apply voltage Vg to an additional back- 
gate capacitively coupled to the QD, see Fig. 3.1(a), the Coulomb blockade can be 
lifted. Indeed, by changing Vg one can change the electrostatics so that energies of 
the quantum dot with n and n + 1 electrons become equal, and so an electron can 
freely jump from the left lead onto the QD and then out to the right lead. Thus, 
a current event has occurred, and a peak in the conductance corresponding to that 
back-gate voltage, V^^„^_i, is observed. By sweeping the back-gate voltage, a series of 
peaks is observed. Fig. 3.1(b). 

In this chapter we calculate the correction to the spacing between Coulomb block- 
ade peaks due to finite dot-lead tunnel coupling. In recent years, low-temperature 
Coulomb blockade experiments have been repeatedly used to probe the energy level 
structure of quantum dots [5, 11, 10]. The dot-lead tunnel coupling discussed here 
may influence such a measurement - the presence of leads may change what one sees 
- and so an understanding of coupling effects is needed. One dramatic consequence is 
the Kondo effect in quantum dots [8, 37, 9]. Here we assume that T':^Tk, where Tk 
is the Kondo temperature, and, therefore, do not consider Kondo physics, focusing 
instead on less dramatic effects that, however, survive to higher temperature. 

We study an ensemble of chaotic ballistic (or chaotic disordered) quantum dots 
with large dimensionless conductance, g ^ 1. The dimensionless conductance is 
deflned as the ratio of the Thouless energy Et to the mean level spacing A: g = 
Et/A [10]. For isolated quantum dots with large dimensionless conductance, the 
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Figure 3.1: (a) Scheme of the Coulomb blockade setup; (b) Oscillations of the 
SET linear conductance G as the back-gate voltage Vg is changed. "Even" ("odd") 
corresponds to an even (odd) number of electrons in the valley. Arrows show average 
shifts in the positions of the peaks' maxima due to the finite dot-lead couplings. 



distribution of g energy levels {sk} near the Fermi level and the corresponding wave 
functions {■^/'^(r)} can be approximated by random matrix theory (RMT) [38, 10, 18]. 
As will be evident from what follows, the leading contribution to the results obtained 
here comes from ^ = 2Ec/A energy levels near the Fermi level; thus, if Ec ^ Et 
the statistics of these ^ levels can be described by RMT. We furthermore neglect the 
spin-orbit interaction and, therefore, consider only the Gaussian orthogonal (GOE) 
and Gaussian unitary (GUE) ensembles of random matrices. 

The microscopic theory of electron-electron interactions in a quantum dot with 
large dimensionless conductance brings about a remarkable result [19, 12]. To leading 
order, the interaction Hamiltonian depends only on the squares of the following two 
operators: (i) the total electron number operator n = Y^ c^ka^ko- where {cka] are the 
electron annihilation operators and a labels spin, and (ii) the total spin operator 
S = I X] '^lo-i {^i\^\^2) Ckcr2 where {ai} are the Pauli matrices. The leading order 
part of the Hamiltonian reads [19, 12] 



H, 



int 



Ecfi^-Js^'' 



(3.1) 



where Eq is the redefined value of the charging energy [19, 7] and Jg > is the exchange 
interaction constant. Higher order corrections are of order /\/ g [19, 12, 20, 7]. The 
coupling constants in (3.1) are invariant with respect to different realizations of the 
quantum dot potential. This "universal" Hamiltonian is also invariant under arbitrary 
rotation of the basis and, therefore, compatible with RMT. In principal, the operator 
of interactions in the Cooper channel can appear in the "universal" Hamiltonian for 
the GOE case. However, if the quantum dot is in the normal state at T = 0, then 
the corresponding coupling constant is positive and is renormalized to a very small 
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value [19, 39]. The "universal" part of the Hamiltonian given by Eq. (3.1) is called 
the constant exchange and interaction (CEI) model [20, 7]. 

The total Hamiltonian of the quantum dot in the g^oo limit thus has two parts, 
the single-particle RMT Hamiltonian and the CEI model describing the interactions. 
The capacitive QD-backgate coupling generates an additional term which is linear in 
the number of electrons: 

i^dot = 5^ ^fe cl^Cka + Ec{n- Mf - Js S^ (3.2) 

kcr 

where J\f=CgVg/e is the dimensionless back-gate voltage and Cg is the QD-backgate 
capacitance. 

The CEI model contains an additional exchange interaction term as compared 
to the conventional constant interaction model (CI model) [40, 41, 5]. Exchange is 
important as Jg is of the same order as A, the mean single-particle level spacing. 
Indeed, in the realistic case of a 2DEG in a GaAs/AlGaAs heterostructure with gas 
parameter Vg = 1.5, the static random phase approximation gives Jg ~ 0.31A [42]. 
Therefore, as we sweep the back-gate voltage, adding electrons to the quantum dot, 
the conventional up-down filling sequence may be violated [21, 22]. Indeed, energy 
level spacings do fluctuate: If for an even number of electrons n in the QD the 
corresponding spacing, eii+i — en, is less than 2Js, then it becomes energetically 
favorable to promote an electron to the next orbital instead of putting it in the same 
one; thus, a triplet state {S = 1) is formed. Higher spin states are possible as well. 
For rs = 1.5 the probability of forming a higher spin ground state is Pi{S>0) ~ 0.26 
and P2{S > 0) ~ 0.19 for the GOE and CUE, respectively. The lower the electron 
density in the QD, the larger r^ and, consequently, the larger the exchange interaction 
constant Jg. 

The back-gate voltage corresponding to the conductance peak maximum A/'„__ i is 
found by equating the energy for n—1 electrons in the dot with that for n electrons:^ 

En-liK^^J = i5n(AC_.) . (3.3) 

As we are interested in the effect of dot-lead coupling on these peak positions, it is 
natural to expand the energies perturbatively in this coupling: E = E^^' + E^"^' + 
.... One possible virtual process contributing to E'^'^'> is shown in Fig. 3.2. Electron 

"'^More precisely, the conductance peak has its maximum at the gate voltage corresponding to 
the maximum of the total amplitude of the tunneling through many energy levels in a QD: A — 
Tlika -^ka, see Ref. [43]. Therefore, in addition to the dominant resonant term which gives the charge 
degeneracy point result for the peak maximum [Eq. (3.3)], the total tunneling amplitude includes 
elastic cotunneling terms (we assume that T <C A) with random phases. However, one can show 
(see Ref. [43] for details) that these cotunneling terms give negligible fluctuating correction to the 
position of the peak maximum with mean((5A/'„_ i) ~ and rms((5A/'„_i) ~ {9l + 9r) / (47r) ^. 
The coefficient in the last equation corresponds to the Cf model. 
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Figure 3.2: One example of the virtual processes contributing to E'^'^\ This virtual 
process corresponds to an electron tunneling out of the quantum dot into the left lead 
and then tunneling back into the same level in the QD. 



occupations of the QD "to the left" and "to the right" of the conductance peak [see 
Fig. 3.1(b)] are different; hence, the corrections E^"^' to the energies are different. 
Therefore, the position of the peak maximum acquires corrections as well. A/" = 

A/''-'^-' + A/"*-^-* + . . . , as does the spacing between two adjacent peaks. 

This physical scenario has been considered by Kaminski and Glazman with the 
interactions treated in the CI model, i.e. neglecting exchange [43]. The ensemble 
averaged change in the spacing and its rms due to mesoscopic fluctuations were cal- 
culated. On average, "even" spacings (that is, spacings corresponding to an even 
number of electrons in the valley) increase, while "odd" spacings decrease (by the 
same amount) [43]: 

2Ec 



/{2) 



2EcU}i'{Js = 0) = A 



9l + 9r 



In 



(3.4) 



27r2 '" T 

where U is the dimensionless spacing normalized by 2Ec and gL^fi = G l,r/ {^e"^ / h) are 
the dimensionless dot-lead conductances. 

In this chapter we calculate the same quantities but with the electron-electron 
interactions in the QD described by the more realistic CEI model. We find that the 
average change in the spacing between conductance peaks is significantly less than 
that predicted by the CI model. However, the fluctuations are of the same order. In 
contrast to the CI result [43], for large enough Jg, we find that "even" spacings do 
not necessarily increase (likewise, "odd" spacings do not necessarily decrease). 

The chapter is organized as follows. In Section 3.3 we write down the total 
Hamiltonian of the system and find the condition for the tunneling Hamiltonian to 
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be considered as a perturbation. In Section 3.4 we describe the approach and make 
symmetry remarks. In Section 3.5 we perform a detailed calculation of the correction 
to the spacing between Coulomb blockade peaks for the I ~^ -'■ ~^ I ^P^'^ sequence. In 
Section 3.6 we find the ensemble averaged correction to the peak spacing. The rms 
of the fluctuations of the correction to the peak spacing is calculated in Section 3.7. 
In Section 3.8 we summarize our findings and discuss their relevance to the available 
experimental data [44, 45]. 



3.3 The Hamiltonian 

The Hamiltonian of the system in Fig. 3.1(a) consists of the QD Hamiltonian 
[Eq. (3.2)], the leads Hamiltonian, and the tunneling Hamiltonian accounting for the 
dot-lead coupling: 

H = Hdot + -f^lcads + Htun- (3.5) 

The leads Hamiltonian can be written as follows 

-"leads = _^ ^p'^pcr'^po- + _^ ^q'^qcr'^qo-) (3-0) 

per qcr 

where {^p} and {^q} are the one-particle energies in the left and right leads, respec- 
tively, measured with respect to the chemical potential (see Fig. 3.2). We assume 
that the leads are large; therefore we (i) neglect electron-electron interactions in the 
leads and (ii) assume a continuum of states in each lead. The tunneling Hamiltonian 
is [46] 

^tun = ^{tkpcl^Cp^ + h.c.) + ^(tfeqC^Cqa + h.c), (3.7) 

fcpcr fequ 

where {t^p} and {tkq} are the tunneling matrix elements. 

We assume that T ^ A and, therefore, neglect excited states of the QD con- 
centrating on ground state properties only. We also assume that the QD is weakly 
coupled to the leads, treating the tunneling Hamiltonian as a perturbation. Correc- 
tions to the position of the peak maximum can be expressed in terms of corrections to 
the ground state energies of the QD via Eq. (3.3). The perturbation series for these 
corrections contains only even powers as i/tun is off-diagonal in the eigenbasis of Hq. 
The 2m*'^ correction to the position of the peak is roughly 



gL + gR, 2EcfgL + gR, 2Ec , 2E, 



c 



m—l 



2Ec Ar(2™) ^ A ^^i^ In ^ ( ^^i^ In ^ In ^ J . (3.8) 

Thus, finite-order perturbation theory is applicable if [43] 

»- + »«,u(?^V.>f^l«l. (3.9) 
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To loosen this restriction one should deploy a renormalization group technique which, 
however, is beyond the scope of this chapter [43, 47, 48, 49]. 



3.4 Plan of the Calculation 

As the exchange interaction constant Jg becomes larger, more values of the QD 
spin S become accessible. The structure of the corrections to the ground state energies 
depends on the total QD spin S, and this structure becomes very complicated for large 
values of the spin S. Fortunately, for the realistic case r^ = 1.5, the probability of 
spin values higher than | in an "odd" valley is small: P2{S> |) ~ 0.01 for the GUE. 
Hence, we can safely assume that in the "odd" valley the spin is always equal to |. 
In the "even" valley, one has to allow both S = and S = l states. 

The structure of the expression for the spacing between peaks depends on the al- 
lowed spin sequences. For an "even" valley there are only two possible spin sequences: 

I^O^I and l^l^l (3.10) 

where the number in the middle is the spin in the "even" valley, while the numbers 
to the left and right are spin values in the adjacent valleys. Fig. 3.1(b). For an "odd" 
valley there are four possibilities: 

0^1^ 0, 0^1^ 1, 
l^i^O, and l^i^l. (3.11) 

To obtain correct expressions for the average spacing between peaks, one should 
weight these sequences with the appropriate probability of occurrence. 

Before proceeding with the calculations, we note several general properties. First, 
ensemble averaged corrections to the "odd" and "even" spacings are of the same 
magnitude and opposite sign. Fig. 3.1(b). Second, the mesoscopic fluctuations of 
both corrections are equal. Indeed, the shift in position of an "even-odd" (n — ^ n+1) 
peak maximum. Fig. 3.1(b), is determined by the interplay between the ^ | and 
1-^1 spin sequences. Likewise, the shift of the "odd-even" (n — 1 -^ n) peak is 
determined by the 5 — ^ and ^ — ^ 1 spin sequences. Now if we sweep the back- 



2 " " ^ 2 
gate voltage in the opposite direction and write the same peak as n — *> n— 1, then the 

corresponding spin sequences are exactly the same as they were in the flrst case: — > | 

and 1 ^ |. From this symmetry argument one can conclude that (i) the ensemble 

averaged shifts of the "even-odd" and "odd-even" peaks are of the same magnitude 

and in the opposite directions [43] and (ii) the mesoscopic fluctuations of both shifts 

are equal. 

Thus, to simplify the calculations we study only the "even" spacing case. This 
corresponds to the two spin sequences given in Eq. (3.10). First, we calculate correc- 
tions to the spacing between peaks for both spin sequences. A complete calculation 
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Figure 3.3: Occupation of tlie QD levels in the ground state in three consecutive 
valleys with total electron number n — 1, n, and n + 1, respectively. A doublet- 
triplet-doublet spin sequence is shown. The variables x, (xi, 0:2, • • • )> ^ind (t/i, 1/2, • • • ) 
denote the energy level spacings in the QD normalized by mean level spacing A. For 
example, x= (eii+i — en) /A. 



for the doublet-triplet-doublet spin sequence is in the next section. Second, we elab- 
orate on how to put these spacings together in the final expression for an "even" 
spacing. Finally, we calculate GOE and GUE ensemble averaged corrections to the 
spacing and the rms fluctuations. 

3.5 Doublet-Triplet-Doublet Spin Sequence: Cal- 
culation of the Spacing Between Peaks 

Let us find the correction to the spacing between peaks for a doublet-triplet- 
doublet spin sequence. The corresponding electron occupation of the quantum dot in 
three consecutive valleys with n — 1, n, and n+1 electrons is shown in Fig. 3.3. 

For the isolated QD the position of the n— 1 -^ n conductance peak maximum is 
determined by 

<°i,.s..(^fi)=<Ui(A^f.) (3.12) 

where E _ 1 and -E„ g^-^ are the ground state energies of dot Hamiltonian [Eq. (3.2)]. 
The corrections due to dot-lead tunneling are different for the doublet and triplet 
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states. The resultant shift in peak position is given by [43] 

•^.f-'i = ^ kk.(A(;°'.) - B;,",,..i(A^f i)l ^ (3.13) 



c 



Note that for the second-order correction to the position, the ground state energies 
are taken at the gate voltage obtained in the zeroth-order calculation, Eq. (3.12). 

Analogous equations hold for the n ^ n+1 conductance peak. The spacing 
between these two conductance peaks is then defined as 

W.,5=i = AC^i (1 - I) - AC_i (I - 1). (3.14) 

3.5.1 Zeroth Order: Isolated Quantum Dot 

For the doublet-triplet n — l^n sequence Eq. (3.12) gives 

Ar(°\ (l^l)=n-i + ^ fe.+i - -J, - - In -^ (3.15) 

"-5^2 ; 2 2Ec V ' 4 2 2/ ^ ' 

where the last temperature dependent term is the entropic correction to the condition 
of equal energies [31]. For the position of the n — * n+1 peak maximum we obtain 

A^fA(l-^)-"4 + ^(^t + J^. + ^'4). (3.16) 

Thus the spacing between peaks in zeroth order is 

«S=i(-) = l + ^ + ^ln|, (3.17) 

where j = Js/A and x = (eii+i — en) /A (see Fig. 3.3). Similarly, for the doublet- 
singlet-doublet spin sequence the spacing is 

W^o) (^) ^ 1 + ^^:i^ - ^ ln2. (3.18) 

n,S-0\ ) 2^ 2Ec 

Note that in both cases Un depends only on the spacing x. 



3.5.2 Second Order: Contribution From Virtual Processes 

Let us consider in detail the second-order correction to the ground state energy 
of the triplet for subsequent use in (3.13): 

<L.(A^) = E l;, ... .1.,/. . (3.19) 



BJS.I W - K. (A/") 
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Figure 3.4: Six distinct types of the virtual processes contribute to the ground state 
energy correction for the QD in the triplet state. Only tunneling processes in (or out 
of) the left lead are shown. In the first four cases spin of the dot in the virtual state S 
has a definite value. In the last two cases (e) and (f ) QD spin in the virtual state has 
two allowed values: S = \ and 5 = 1 with the probabilities ws given by Eq. (3.21). 
Electron structure of the virtual state corresponding to two allowed values of S is 
circled by dashed line in panel (e). 
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where the sum is over all possible virtual states. E^'^^ and |-?/''^°^) are the eigenvalues 
and eigenvectors of H^ot, Eq- (3.2). 

Different terms in Eq. (3.19) have different structure depending on the type of 
virtual state involved; six possibilities are shown in Fig. 3.4. To take into account all 
virtual processes, we sum over all energy levels in the QD and integrate over states 
in each lead. To simplify the calculation even further, we assume (just for a moment) 
that T = so that the Fermi distribution in the leads is a step function. Later we will 
see how T reappears as a lower cutoff within a logarithm. 

Following the order of terms in Fig. 3.4, the second-order correction to the triplet 
ground state energy is 

oo Of \ \+ |2 



E^S-m = - E E 



Ek - Sp) + 2Ec{n + 1 - AT) - I J3 



h2 p ^ " P^ — ov- 1 2 ■" ^ 4 
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'2~^ nl \ u |2 



(^p) IWl 



EE 

fe = l p V-P -K7 — OV" 2 ■'-I 4 



ep-ek)-2Ec{n-\-M)-\J, 



r+l 



- EE 

fc=f p 
f+i 

- EE 



l^^p/ |''A;p| 



,= . p ,ep-5fc)-2Ec(n-i-Ar) + |j, 

f + 1 /)/ N I , |2 

f^l^^pj Kfcpl 
,= . p (£fc-ep)+2Ec(n+|-Ar) + fJ3 






tt V (^P - ^^) - 2^^(^ - 2 - ^) + /^-^^ 



|2 

+ E 



■^^ 6'(— ep) |tfcp| 

,=:,+2 p (^it-^p)+2i?c(n+|-AA) + /5Js 



Z^ 



2 



+ {similar terms for the right lead: p ^ q} , (3.20) 

where S is the spin of the QD in the virtual state. One can easily find S for the first 
four processes. Figs. 3.4(a)-(d), and so calculate the denominators for the first four 
terms in (3.20): the values are |, |, |, and |, respectively. In the last two cases. 
Figs. 3.4(e)-(f), the QD spin in the virtual state can take two values, 5 = | or |; it 
does so with the following probabilities 

wi = I and wz = \. (3-21) 

The corresponding contributions to the energy correction must be weighted accord- 
ingly. In addition, the energy difference in the denominators depends on S] to account 
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for this dependence, we introduce an additional function 

/5 = 2-5(5 + l), 



(3.22) 



appearing in the denominators of the fifth and sixth terms in (3.20). 

Let us integrate over the continuous energy levels in the lead. The sum can be 
replaced by an integral, ^ ■ ■ ■ — > f de ul^e) ■ ■ ■ , where z/^ is the density of states 
in the left lead. Taking the dot-lead contacts to be point-like, the tunneling matrix 
elements {tkp} depend on the momentum in the lead p only weakly; hence, tkp ~ tki- 
In addition, as the leads are formed from 2DEG, their density of states is roughly 
independent of energy. We assume that it is constant in the energy band of 2Ec near 
the Fermi surface. Then the result of integrating over the energy spectrum in the 
lead (in schematic form) for the first term in Eq. (3.20) is 



p 






J^L \t 



kL\ 



In 



(3.23) 



This expression diverges, but when we calculate an observable, e.g. the shift in the po- 
sition of the peak maximum [Eq. (3.13)], we encounter the energy difference between 
corrections to the triplet and doublet states. The result for the shift is, therefore, 
finite: 
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Cfe 



(3.24) 



In a similar fashion one can calculate the second-order correction to the ground 
state energy of the doublet. The difference of these energies at the gate voltage 
corresponding to the peak maximum in zeroth order [Eq. (3.15)], needed in Eq. (3.13), 
then follows. There is one resonant term, proportional to ln(2i?c/T), in which the 
lower cutoff T appears because of the entropic term in Eq. (3.15). Alternatively, T 
would appear as the natural cutoff for the resonant term upon reintroduction of the 
Fermi-Dirac distribution for the occupation numbers in the leads. 

For a point-like dot-lead contact, the tunneling matrix element is proportional to 
the value of the electron wave function in the QD at the point of contact: t^^ oc ipki^a)) 
where a = L or R. Here, we neglect the fluctuations of the electron wave function in 
the large lead. Thus, the following identity is valid 



^a l^-fec 



A \iJk{ra)\ 

9a 



47r2 



2\ ' 



[\^kira)\ > 



(3.25) 



where the average in the denominator is taken over the statistical ensemble. Note that 
by taking the ensemble average of both sides of (3.25), one arrives at the standard 
golden rule expression for the dimensionless conductance. 
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In our calculations we take advantage of the fact that Jg < A ^ Ec and neglect 
terms that are of order l/^ = A/2Ec. Sums like 






2X 



Ek - £§+1 



(3.26) 



r+?+l 



are split using Xlfel^+2 ' ' ■ = X]fc=ii+2 ' ' ■+X]fcl2+f+2 ' ' ' ' ^^"^ ^*-' ^^^ ^^^^ term, which 
is 0(1/^), is dropped. Likewise, expressions like 



2 ^ , / 3 



J. 



fc=f+2 



2 Ek-ei 



+1 



2E, 



c 



(3.27) 



are of order 0(1/.^) and so neglected. 

Thus, for the second order correction to the position of the peak maximum, we 
obtain 
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(3.28) 



where 2 Jg > ^n+i— e^ > because the total spin of the QD with n electrons is equal 
to 1. 

In a similar fashion one can find the shift in the position of the other peak 
maximum, A/li^ {^ ^ ^)- Then, according to (3.14), the difference of these two 
shifts yields the second-order correction to the spacing for the doublet-triplet-doublet 



48 



spin sequence: 
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(3.29) 



A potential complication is that the addition of two electrons to the quantum dot 
(n — 1 — ^ n -^ n + 1) may scramble the energy levels and wave functions of the 
QD [50, 51, 7]. Since the number of added electrons is small, we assume that the 
same realization of the Hamiltonian, Eq. (3.5), is valid in all three valleys [52, 53]. 

For the second-order correction to the spacing for the doublet-singlet-doublet spin 
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sequence we similarly obtain 
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(3.30) 



where snj^i—sn. > 2 Jg > because the total spin of the QD with n electrons is equal 
to in this case. 

Unlike the zeroth-order spacings, the second-order corrections are functions of 
many energy level spacings as well as the wave functions at the dot-lead contact 
points: W;^^ = U:^'s{x,X,Y; {Zka}), where x, X = {xi,X2,...), and Y = {yi,y2,---) 
are the energy level spacings in the QD normalized by the mean level spacing A (see 
Fig. 3.3) and 

l^fc(ra)|^ 



Zhn = 



[\^k{ra)f) 



(3.31) 



The expressions for W*^^-* suggest that the main contribution to their fluctuation 
comes from the fluctuation of the energy level x and the wave functions {ipk (i"a)}. 
The other spacings, X and Y, always appear within a logarithm; therefore, their 
contribution to the fluctuation of W*^^-* is small. With good accuracy, one can replace 
these levels by their mean value 

wi5(x,l,l;{Z,J)=wS(x;{Z,J). 



n,S 



u'-' - 



Converting to dimensionless units, we flnd that 
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(3.32) 



(3.33) 
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where 



^a,S=0 {x; {Zka}) 
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(3.35) 



where (5 = A/T. Here, the upper hmit in two of the sums is infinity because the 
Fermi energy is the largest energy scale. 

In summary, the total spacing is 



/(O) 



(2), 



U^,s = K;s{x) +Ui:;s(.x-{Zka}) 



(3.36) 



where the first term is given by Eqs. (3.17) and (3.18) and the second by (3.33)-(3.35). 
The spin of the QD in the even valley, S, can take two values, or 1, depending on 
the spacing x. 



3.6 Ensemble Averaged Correction to the Peak 
Spacing 

The average and rms correction to the peak spacing can now be found by using 
the known distribution of the single-particle quantities x and {Z^a}- In what follows, 
(U) denotes the average over the wave functions, U denotes the full average over both 
wave functions and energy levels, and P{x) is the distribution of the spacing x. Since 
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(Zfca) = l, {^a,s) does not depend on a, and the average "even" spacing is^ 



Ui'^ = j^ dxP{x)(u^,) (3.37) 



^^ (^j^ dx P{x) ($5^0) + j^x P{x) ($S=i)) • (3.38) 



Using the asymptotic formulas 






X In^, 



E ['"(^ + t)-'"U^ + ,;j 

i = l 

for ,^^1 in the expressions for ($^5), we find 

($«,s=o) = 2(2a;-3j-l)lne + 2 1n(5, 

{^a,s=i) = 2{-2x + 5j-l)\n^-\n6, (3.40) 

vahd for ^,S ^ 1. By carrying out the integration over the distribution of the spacing 
X, the final expression is 

l(F(j) = ^-^^[Cij)\n^ + Vij)\n5 + 0(1)1 

C{j) = 2[-8jPo{2j)+4xo{2j) + 5j-3], (3.41) 

Vij) = 3Po{2j)-l, 

where Po{2j) = f^ dx P{x) and xq{2j) = f^ dxx P{x). Note that Po{2j) is the 
probability of obtaining a singlet ground state while Xq{2j) / Pq{2j) is the average 
value of X given that the ground state is a singlet. 

For the CI model, j = and, hence, C(0) = 1^(0) = 2. In this limit, then, the 
ensemble averaged correction to the spacing is 

77(2)777 9l + gR , 2Ec , . 

in agreement with previous work [43]. The magnitude here is approximately 0.05 
{gi + gn) In {2Ec/T) in units of the mean level spacing. 



^One may think that an additional term due to shift in the border x = 2j between singlet and 
triplet states should be present in this equation as well. However, one can show that the contribution 

(9) 

of this term to Un (as well as to its rms) is small. 
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Figure 3.5: Correction to "even" peak spacing as a function of strength of exchange, 
j = Js/A, normalized by the correction at j = 0. GUE case with A = 1 and gi = Qr- 
Sohd: Ensemble average. Dashed: Ensemble average plus/minus the rms, showing 
the width of the distribution. 



Js 



It is convenient to relate the average change in spacing at non-zero J^ to that at 
= 0: 



5u{j) 
A 



U^n^ (j) 



XC{j)+V{j) 



Ine 
In (5 



(2) t^^ 2(A + 1) 

ln(2Ec/A) 



W (0) 



ln(A/r) 



(3.43) 
(3.44) 



The dependence of 5u on j is fully determined by the parameter A and the choice of 
random matrix ensemble. 

Figs. 3.5 and 3.6 show the results in the GUE ensemble for A = l and 3, respec- 
tively, and Fig. 3.7 shows those for the GOE ensemble at A = 1. In evaluating these 
expressions, we use the Wigner surmise distributions for P{x), which allow an analytic 
evaluation of Po{2j) and xo(2j). As j increases, the average correction to the peak 
spacing decreases monotonically in all three cases. (Note, however, that our results 
are not completely trustworthy when 0.4<j<0.5 because in this regime higher spin 
states should be taken into account.) Since A depends on ^ and 6 only logarithmically, 
the qualitative behavior of Su{j) is very robust with respect to changes in charging 
energy, mean level spacing, or temperature. 
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Figure 3.6: The same quantities as in Fig. 3.5 plotted for A = 3. 
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Figure 3.7: Ensemble averaged correction to the "even" peak spacing as a function 
of strength of exchange, j = Jg/A, normalized by the correction at j = for A = 1. 
Solid: GOE. Dotted: GUE. 
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Figure 3.8: GUE ensemble averaged correction to tlie even peak spacing as a function 

of A = ln(2£'c'/^)/ln(A/T) at j=jo- The curves correspond to jo = 0.3, 0.4, and 0.5 
from top to the bottom. 



Similarly the dependence of 6u on A at j =jo is 

\C{jo)+Vijo) 



Su{jo, A) 



2(A + 1) 



(3.45) 



Fig. 3.8 shows results in the GUE case for several values of jo- Thus for the realistic 
value jo = 0.3, the CEI model gives an average correction to the peak spacing which is 
two to three times smaller than the CI model. 



3.7 RMS of The Correction to Peak Spacing due 
to Mesoscopic Fluctuations 

Mesoscopic fluctuations of the correction to the peak spacing are characterized by 
the variance oiU^'^\ It is convenient to separate the average over the wave functions 
from that over the spacing x, writing 



where 



var(wf)=4(Wf)+a^(Wf), 

poo 

4 = I c/a:P(x)((wS-(wS> 



(3.46) 
(3.47) 
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is the contribution due to wave function fluctuations and 



/■oo 

al= / rfxP(x)(wS>^- 
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dxP{x){U)^^l) 



-[ 2 
(2)' 



(3.48) 



is the contribution due to fluctuation of the spacing x. 

We start by considering the fluctuations of the wave functions. As the number of 
electrons in the dot is large, the distance between the left and right dot-lead contacts 
is large, jr^ — r^jl ^ \p where Xp is the Fermi wavelength. Therefore, the wave 
functions at yl and yr are uncorrected [54], 

{{ZuL - 1) {Zk'R - 1)) = (3.49) 

for all k and k' . The fluctuation of W*-^-* can then be written entirely in terms of the 
properties of a single lead: 

((«S - («S>)') = |^((*-^- - <*--) f) (3.50) 

The cross terms here disappear because, according to RMT, wave functions of different 
energy levels are uncorrelated even at the same point in space, 

{{ZkL-l){Zk'L-l)) = '^kk' (3.51) 

where 13 = 1 {13 = 2) for the GOE (GUE) case. In fact, only the A; = f and A; = f + 1 
terms contribute, as one can see by using 



5:'"1i-^)-«w f^'^^) 



valid for .^ ^ 1. Integrating (3.50) over the distribution of x according to Eq. (3.47) 
(keeping in mind ^,5^1), we obtain 

^^(^'^0 = ^TT^i^l^'^ + [3^o(2j) + 1] In^^ - 4 [3Po(2j) - 1] \H\n5]. (3.53) 
^ ^ [3 (47r^^) 

In the contribution to the variance due to fluctuation of the level spacing x, 
Eq. (3.48), (W„ g) can be taken from the previous section. Since the average eliminates 
the dependence on the lead a, we have immediately 

{l^Z'sf = [^)\^^.f (3.54) 
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where {^l,s=o) and (<l>i 5=1) for ,^,5^1 are given by Eq. (3.40). Using these expres- 
sions in Eq. (3.48), we obtain 

al = (^^^) (Qs In' e + Css In^ 5 + C^s In^ In 5) . (3.55) 

Exphcit expressions for the coefficients are given below once we reach the final result. 

The dependence on gi and qr of the two contributions to the variance is different. 
In particular, the contribution due to fluctuations of the wave functions [Eq. (3.53)] 
is proportional to 

<,i + jj=(££±i^+fe£^, (3,56) 

The first term has the same form as the contribution (3.55) from fluctuations of x. 
It is convenient to write the total variance as a sum of symmetric and antisymmetric 
parts. Our final result for the variance is 

var(Wf)=af(Wf)+a^(Wf) (3.57) 

where 

^.'(Wf) = (^^^^J[s,,ln'^ + Sssln'6 + S,sln^ln6) (3.58) 

^'(Wf) = (^1^) {Aiiri'^ + Assln'6 + A^sln^\n6) (3.59) 



with the coefficients {S} and {A} given by 

%(j) = | + 16(x-l) 

+ 64 [2jPo(2j) - xo{2j)] {2j [1 - Po(2j)] " [1 " :^o(2j)]} , (3.60) 

SssU) = ^[3Po(2j) + l] + 9Po(2j)[l-^o(2j)], (3.61) 

Mj) = -|[3Po(2j)-l] + 24{xo(2j)[l-Po(2j)] 

+ [l-xo(2j)]Po(2j)-4jPo(2j)[l-n(2j)]}, (3.62) 

Adj) = |, Assij) = ^[SPo{2j) + l], ^e,(j) = -|[3Po(2j)-l]-(3.63) 

The constant x introduced in Eq. (3.60) is 

POO 

X= dxx^P{x). (3.64) 

Jo 
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Figure 3.9: The rms of the correction to the peak spacing for the symmetric setup 
Ol = Or diS a. function of j = Jg/A normahzed by the ensemble averaged correction 
at j = 0, Eq. (3.66). Sohd (dotted) curve corresponds to the GOE (GUE) at A = 1. 
Dashed curve corresponds to the GUE at A = 3. 



-'("^')U=|7fe|('.>«-ln^)Vl6(x--l)(^) ln=«. (3.65) 



For the CI model, j = and Po(0) = Xo{0) = 1; hence 

il + g' ' ^' 

a (47r20 

The first term is due to fluctuation of the wave functions at the dot-lead contacts, 
and the second term comes from the fluctuation of the level spacing x. The presence 
of the second term was missed in previous work [see Eq. (44b) in Ref. [43]]. If .^ = (5, 
the first term vanishes; nonetheless, due to the second term the variance is always 
finite. 

Let us consider a realistic special case of symmetric tunnel barriers, gi, = qr [44]. 
Then the asymmetric contribution vanishes, and the rms fluctuation of the correction 
to the peak spacing normalized by the average correction at j = [Eq. (3.42)] is 

. (A - ^-(^"'^) - VMJ)>^' + Mj)>^ + Sssij) ._.. 

Ui^^ (0) 2(A + 1) 

Fig. 3.9 shows this quantity plotted as a function of j for both GOE and GUE. Notice 
that (i) the rms is of the same order as the average, and (ii) its magnitude weakly 
depends on j. To show the magnitude of the fluctuations in the correction relative to 
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its average value, we plot two additional curves in both Fig. 3.5 and Fig. 3.6, namely 
5u ± asu- We find that at the realistic value j = 0.3, the correction to the even peak 
spacing is negative for a small fraction of the quantum dots in the ensemble. 

3.8 Conclusions 

In this chapter we studied corrections to the spacings between Coulomb blockade 
conductance peaks due to finite dot-lead tunneling couplings. We considered both 
GOE and GUE random matrix ensembles of 2D quantum dots with the electron- 
electron interactions being described by the CEI model. We assumed T^A^i^^- 
The S* = 0, |, and 1 spin states of the QD were accounted for, thus limiting the 
applicability of our results to Js<0.5A. 

The ensemble averaged correction in even valleys is given in Eq. (3.41). The 
average correction decreases monotonically (always staying positive, however) as the 
exchange interaction constant Jg increases (Figs. 3.5-3.7). The behavior found is very 
robust with respect to the choice of RMT ensemble or change in charging energy, 
mean level spacing, or temperature. Our results obtained in second-order perturba- 
tion theory in the tunneling Hamiltonian are somewhat similar to the zeroth-order 
results [20, 7] in that the exchange interaction reduces even-odd asymmetry of the 
spacings between peaks. While the average correction to the even spacing is positive, 
that to the odd peak spacing is negative and of equal magnitude. 

The fluctuations of the correction to the spacing between Coulomb blockade peaks 
mainly come from the mesoscopic fluctuations of the wave functions and energy level 
spacing x in the QD. The rms fluctuation of this correction is given by Eqs. (3.57)- 
(3.64). It is of the same order as the average value of the correction (Figs. 3.5 and 
3.6) and weakly depends on Jg (Fig. 3.9). Therefore, for a small subset of ensemble 
realizations, the correction to the peak spacing at the realistic value of j = 0.3 is of 
the opposite sign. The rms fluctuation of the correction for an odd valley is the same 
as that for an even one. 

We are aware of two experiments directly relevant to the results here. First, in the 
experiment by Chang and co-workers [45], the corrections to the even and odd peak 
spacings due to finite dot-lead tunnel couplings were measured. It was found that the 
even (odd) peak spacing increases (decreases) as the tunnel couplings are increased. 
This is in qualitative agreement with the theory, see Eq. (3.41). The magnitude of 
the effect was measured at different values of the gas parameter r^ (and, hence J^) 
as well. Unfortunately, because the effect is small and the experimentalists did not 
focus on this issue, one can not from this work draw a quantitative conclusion about 
the behavior of the correction to the peak spacing as a function of Jg. 

Second, in the experiment by Maurer and co-workers [44], the fluctuations in the 
spacing between Coulomb blockade peaks were measured as a function of the dot-lead 
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couplings with gL = gR- Therefore, only the symmetric part [Eq. (3.58)] would con- 
tribute to the total variance. Ref. [44] reported results for two dots: a small one with 
area 0.3 ;um^ and a large one with area 1/im^. From the area (excluding a depletion 
width of about 70 nm), we estimate that the large (small) dot contains about 500 (100) 
electrons. Measurements on the large QD found larger fluctuations upon increasing 
the dot-lead tunnel coupling, in qualitative agreement with the theory. Though the 
temperature was larger than the mean level spacing in the large QD whereas our 
theory is developed for T<^ A, the theory gives about the correct magnitude for the 
peak spacing fluctuations. It is inconclusive whether the data is in better agreement 
with the CI or CEI model as the fluctuations are roughly the same (Fig. 3.9) in both. 
In the small QD, there is an anomaly for the strongest coupling in the experiment 
- the fluctuations suddenly decrease. In addition, the experimental fluctuations are 
one order of magnitude larger than the theoretical estimate [Eqs. (3.66) and (3.42)]. 
The reason for this discrepancy is not clear at this time. Possible contributing factors 
include scrambling of the electron spectrum as the charge state of the dot changes, 
or the role of the fluctuations when the dot is isolated (i.e., fluctuations in W*^°^). In 
order to assess quantitatively the role of dot-lead coupling in the Coulomb blockade, 
further experiments are needed. 



Chapter 4 

Introduction to Quantum 
Computation with Quantum Dots 

4.1 Overview 

This chapter serves as an introduction to the field of quantum computation (QC). 
The plan of the chapter is as follows. In Section 4.2 we explain why QC field deserves 
our attention. In Section 4.3 we outline the mathematical basis of quantum compu- 
tations. In Section 4.4 we list and briefiy explain DiVincenzo's criteria of physical 
implementation. Finally, in Section 4.5 we briefly mention liquid state NMR quantum 
computing proposal but mostly concentrate on the spintronic quantum dot proposal 
put forth by Loss and DiVincenzo [55]. 

4.2 Introduction to Quantum Computation 

Quantum computation (QC) holds out tremendous promise for efficiently solv- 
ing some of the most difficult problems in computational science: integer factoriza- 
tion [56], discrete logarithms, and modeling of quantum mechanical systems [57]. 
These problems are intractable on any present or future conventional computer. 

The fact that quantum systems can perform a computation was realized in 1982 
by Paul Benioff and Richard Feynman. In 1994 Peter Shor proved that quantum 
computer can perform prime factoring in polynomial time as compared to the classical 
one which requires exponential time to solve this problem [56]. The prime factoring 
problem is defined as follows: given number N which we know is a product of two 
primes P and Q, find P and Q. Table 4.2 gives one a fiavor on how computation 
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Quantum Computer 

of 70 Qubits 
(Shor's Algorithm) 


Classical Computer: 

Network of 

Hundreds Workstations 


130-Digits Number 
Factoring 


1 Month 


1 Month 


400-Digits Number 
Factoring 


3 Years 


10^° Years 
(Age of The Universe) 



Table 4.1: Quantum computer vs. network of hundreds classical computers: Scaling 
of the computational time with the complexity of the factoring problem. 



times scale with the complexity of the factoring problem for classical and quantum 
computers. As one can see for the 400-digits prime factoring problem even relatively 
small quantum computer will do much better than a network of hundreds classical 
computers. This is attained due to quantum parallelism discussed in the next section. 

Shor's spectacular result attracted to the field of QC many physicists and com- 
puter scientists alike. Many proposals for physical realization of QC have been put 
forth since then. New experimental techniques were developed to control phase co- 
herence of the quantum mechanical systems on a nano scale. 

Since prime factoring is the most secure encryption scheme in wide use today, the 
field of QC enjoyed significant funding from government agencies in the USA, EU, 
and Austraha. 



4.3 Mathematical Basis of Quantum Computations 

In this section we review the mathematical basis of quantum computations. 



4.3.1 Quantum Bits and Entangled States 

The bits of information in a quantum computer are formed by two orthogonal 
states of a two-level quantum mechanical system |0) and |1), Fig. 4.1. To distinguish 
between classical and quantum bits quantum bit was given a special name: qubit. For 
example, a nucleus of spin-| (or an electron) can serve as a natural qubit: |0) = |i), 
|1) = It). In contrast to the classical bit, qubit has the following properties [58]: 
(i) The states |0) and |1) do not constitute all possible states of a qubit. In fact, it 
can be in an arbitrary superposition state: |x) = a\0) + (3 |1), where a and P are 
complex numbers with |a| + \f3\ =1 constraint. 
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Figure 4.1: Schematic picture of tlie quantum bit. 



(ii) Entangled state of two qubits is possible, e.g.: 

|z^) = i=(|0)i®|0)2 + |l)i®|l)2). (4.1) 

In this state each individual qubit fails to have any state of its own. The information 
is, therefore, encoded only in the correlation between the qubits 1 and 2. Two qubits 
in an entangled state are also called EPR pair [59]. 



4.3.2 Steps in Quantum Computation and Quantum Gates 

To perform quantum computation we should [58] 
(i) Prepare a collection of n qubits in the appropriate initial state at t = 0: \ip{0)). 
The initial state is usually taken to be 

|V^(0)) = |0),®|0)2® ...®|0)„ (4.2) 

or a state constructed from this one by application of the appropriate unitary trans- 
formations. 

(ii) Let that state evolve under the appropriate unitary transformation U. U can 
be realized as the ordinary time evolution of a physical system with appropriately 
constructed time-dependent Hamiltonian H{t): 

m)) = U mO)) = Te-^r.'^rHir) |^(Q)^ _ (4^3) 

(iii) Extract the result of the computation by appropriate measurement on the final 
state. Typical final measurement is made on a subset of individual qubits. 

Among all possible unitary transformations there are those acting on a single 
qubit, a pair of qubits, a set of three qubits, etc. Let us call unitary transformations 
acting only on a single qubit or a pair of qubits quantum gates. In 1995 Barenco and 
coworkers proved that any quantum computation algorithm can be decomposed in a 
series of quantum gates: single qubit and two-qubit XOR gates [60] . 
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4.3.3 Quantum Parallelism 

The computational basis of n-qubit quantum computer is formed by 2" states 
like this one: 

|0)i ® 10)2 ® 10)3 ® 10)4 ® . . . ® |0)„_3 ® |1)„_2 ® Il)n-1 ® Il)n ' (4-4) 

Each qubit in this ket-vector can be in one of the two possible eigenstates. Let 
us denote ket-vector like the one in Eq. (4.4) simply as \x), where x is the integer 
between and 2" — 1, whose binary expansion gives the corresponding ket-vector. For 
example, in the case of Eq. (4.4) the binary expansion is 0000 . . . 0111. It corresponds 

to \x) = |7). 

The concept of quantum parallelism can be illustrated as follows [58] . Consider 
certain two-qubit entangled state and act on it with a quantum gate U. Since unitary 
transformations are linear we obtain: 

u (^ \o), ® 10)2 + ^J^\l), ® |i) J = -^ w |0), ® 10)2 + yl w |i), ® 11)2 . (4.5) 

Two important remarks should be made about Eq. (4.5): (i) as a result of quan- 
tum computation we obtained the same superposition in the final state as we had in 
the initial state and (ii) we performed quantum computation on all input states si- 
multaneously! This property of quantum computation is called quantum parallelism. 
Quantum parallelism gives quantum computer the potential to do tricks that a clas- 
sical computer can perform only with significantly greater computational effort. 

Suppose we have a unitary transformation Uf, associated with a function /. 
Function / replaces a computational basis state |a;) by |/(a;)), where f{x) can take 
the following integer values: f{x) = 0, 1, 2, . . . , 2" — 1. Since Uf is a unitary trans- 
formation, it is invertible. Hence, function / is bijective - its action is nothing but a 
permutation of integers {x}. This restriction on the nature of / is excessively strong. 
Therefore, let us double the number of qubits and define Uf on the 2n-qubit space as 
follows: 

\x) ® |0) ^ \x) ® \f{x)) . (4.6) 

Let us call the first set of n-qubits input register and the second set of n-qubits output 
register [58]. 

To illustrate the concept of quantum parallelism even further let us take the input 
register to be in the state 

2" — 1 

1^) = ^ E 1^) = 4 (I0)i + Il)i) ® • • • ® 4f (I0)n + |1)J • (4-7) 

x=0 * * 
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This input state can be obtained from the state in Eq. (4.2) by subjecting each qubit 
to the one-qubit gate that "rotates" each "spin" 90 degrees about the t/-axis. Under 
the action of Uf 

1^) ® |0) = ^ E 1^) ® |0) ^ ^ E 1^) ® l/(^)) • (4.8) 

x=0 z=0 

Notice what just happened: in one run we calculated f[x) on all 2" possible inputs. 
If we would have 100 qubits in our input register, then f{x) would have been simulta- 
neously calculated on 2^"^° ^ 10^° different inputs. This massive quantum parallelism, 
is well outside the reach of any classical computer. 



4.3.4 Important Remarks 

Useful information has to be extracted from the output of the quantum compu- 
tation, see Eq. (4.8). To illustrate how the measurement process works let us perform 
a measurement on the final state in Eq. (4.8). Namely, let us measure the individual 
states of all qubits in the input and output registers. With the equal probability of 
2~" one can find any state of the input register. Let us assume that we find |xo) as a 
result of the measurement on the input register. Then the result of the measurement 
on the output register is |/(xo)). Thus, we calculated the value of /(xq) at the ran- 
domly chosen xq. All other information has been lost because we collapsed the whole 
superposition in Eq. (4.8) onto a single component |a;o) ® |/(a;o)). 

If Eq. (4.8) was the only quantum algorithm one could come up with, one would 
be no better off than just doing single classical computation on a randomly chosen 
input. Therefore, the goal of designing quantum algorithms is to have easily accessible 
data of interest in the final state and to make sure that the data can not be generated 
by a classical algorithm without enormously greater computational effort. The first 
algorithm which satisfies these conditions was created by Peter Shor for the prime 
factoring problem. 

The other difficulty is to maintain the coherence of an entangled state of many 
qubits. Unfortunately, it is impossible to completely isolate our qubits from the 
environment during a quantum computation. The coupling to the environment leads 
to decoherence or collapse of the quantum state which is especially severe for the 
entangled states of many qubits. The solution to this challenge is to use quantum 
error-correcting codes [61]. We know that these codes can work only for the sufficiently 
small decoherence rates. Therefore, one of the main problems in physical realization 
of the quantum computer is how to reduce decoherence due to different couplings to 
the environment. 
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4.4 DiVincenzo's Criteria of Physical Implemen- 
tation 

Five DiVincenzo's criteria [62] serve as a general guidance to the physical imple- 
mentation of a working scalable quantum computer. They naturally follow from the 
mathematical basis of QC discussed in the previous section and can be formulated as 
follows: 

(i) A collection of the quantum mechanical two-level systems (qubits) is needed. 
Each qubit should be separately identifiable and externally addressable. One should 
be able to add qubits at will. 

(ii) It should be possible to completely (or with very high accuracy) decouple qubits 
from one another. One should be able to set the state of each qubit to |0) in the 
beginning of each computation. 

(iii) Long decoherence times are required. For the error correcting-codes to be ap- 
plicable decoherence time should be at least 4 orders of magnitude larger than the 
"clock time" of the quantum gates. 

(iv) Logic operations should be doable. For example, two-body Hamiltonians involv- 
ing nearby qubits should be under independent and precise external control. One 
should be able to smoothly (on the time scale of one clock cycle) turn these Hamil- 
tonians on and off. The integral of the pulse should be controlled with the accuracy 
of at least 1 part in 10'^. 

(v) Projective quantum measurements on the qubits must be doable. It is useful, 
though not necessary, for these measurements to be doable fast (that is, within few 
clock cycles) and have high quantum efficiency (otherwise quantum computation has 
to be done in an ensemble style). 

4.5 Quantum Computing Proposals 

4.5.1 Liquid State NMR 

Since the original work by Peter Shor many quantum computing proposals were 
put forth in different areas of physics. At the time of writing liquid state nuclear 
magnetic resonance (NMR) proposal is well ahead of its competition. In the exper- 
iment by Vandersypen and coworkers [63] the custom-synthesized molecule (used as 
the quantum computer) contained five ^^F and two ^^C spin-| nuclei as qubits. The 
quantum gates were realized by a sequence of spin-selective radio-frequency pulses. 
The readout was performed via NMR spectroscopy. For the first time the simplest 
instance of Shor's algorithm: factorization of 15 (15 = 3x5) was demonstrated [63]. 

Unfortunately, liquid state NMR proposal lacks scalability, see the first DiVin- 
cenzo's criterion. Indeed, to realize more complicated algorithms one has to design 
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■ back gates 



magnetized or 
high-g layer 



heterostructure 
quantum well 



Figure 4.2: Schematic picture of tlie spintronic QD proposal by Loss and DiVin- 
cenzo. 2D electron gas is confined on the interface of AlGaAs/GaAs heterostructure. 
The electrodes deposited on top of the heterostructure control the shape of the QDs 
and their tunnel couplings to each other. Back gates control the electron numbers in 
each dot. Magnetized layer is needed to realize one-qubit operations. Copyright by 
G. Burkard, H.-A. Engel, and D. Loss, from Ref. [67]. 



and synthesize new molecules. 



4.5.2 Spintronic Quantum Dot Proposal by Loss and DiVin- 
cenzo 

The experimental progress in solid state QC proposals falls significantly behind 
liquid state NMR proposal. Nonetheless, once the universal gates (one- and two-qubit 
quantum gates) become sufficiently error free the scaling of the quantum computer 
is going to be more or less straightforward. Among different solid state QC pro- 
posals one can distinguish three major categories: (i) superconducting microcircuits 
containing Josephson junctions [64, 65], (ii) nuclear spins of the impurity atoms in 
semiconductor QDs [66], and (iii) electron spin- or charge-qubits in semiconductor 
QDs, see Chapters 5 and 7 for the references. Charge qubit in the semiconductor 
double QD are discussed in Chapter 7. 

Here we describe the original spintronic quantum dot proposal by Loss and Di- 
Vincenzo [55, 62]. They consider an array of coupled 2D lateral quantum dots, see 
Fig. 4.2. Each QD should contain odd number of electrons so that spin of an electron 
on the last occupied orbital serves as a natural qubit: |0) = ||) and |1) = If). 

Let us briefly go over five DiVincenzo's criteria to see if this system satisfies 
them and what is needed to be improved as compared to the modern experimental 
technologies. 

Criterion (i). We have a collection of two-level quantum mechanical systems 
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and we can add to the system as many new qubits as we want. 

For the qubit to be separately identifiable the corresponding quantum dot should 
contain odd number of electrons with sufficiently high accuracy. This electron number 
can be adjusted by back-gate electrode voltages shown in Fig. 4.2. To make sure that 
the number of electrons is odd and hence spin of the QD is equal to | one can employ 
zero-bias anomaly measurement in transport through the QD which signals a Kondo 
effect. 

Criterion (ii). To set the state of each qubit to |0) it is sufficient to place the 
spins in the magnetic field of several Tesla at the liquid-He temperature. For this 
|0) ® |0) ■ ■ ■ |0) state to be stable upon removal of the magnetic field the qubits 
must be decoupled from one another to a very high accuracy. This can be achieved 
by setting the pairs of electrodes between neighboring QDs to high negative voltages 
so that the orbital wave-functions overlap is negligible. 

Criterion (iii). The experiments on the precession of the electron spin in a 
variety of semiconductor materials were performed by Kikkawa and Awschalom [68]. 
In some structures they found spin decoherence times r^ up to hundreds of nanosec- 
onds, that is, at least one order of magnitude larger than charge decoherence times. 
These results are quite encouraging. However spin decoherence times are expected to 
be very device- and structure-specific. Therefore, more experiments especially on the 
QD array setups are necessary. 

One must emphasize that what really matters is the ratio t^Jt^^ where r^ denotes 
the clock cycle time which will be defined below. 

Criterion (iv). To physically realize one-qubit quantum gates it must be pos- 
sible to subject a specified qubit to the localized magnetic field of proper direction 
and strength. This can be done by using a scanned magnetic particle or by the use of 
magnetized barrier material that the electron can be inserted in and out of by electric 
gating, see Fig. 4.2. The experimental realization of both of these ideas is going to be 
quite challenging, especially, taking into account the fact that the required accuracy 
is 1 part in 10*^. 

Fortunately, according to the paper by DiVincenzo and coworkers [69] one can 
completely eliminate one-qubit quantum gates at the expense of increasing the size 
of the logical qubit from one to three QDs. In this scheme two-body Hamiltonians is 
the only building block necessary to realize quantum computations. 

These two-body Hamiltonians, or two-qubit quantum gates if we go back to the 
original Loss and DiVincenzo proposal, can be realized by increasing the voltages 
on the electrodes between neighboring QDs. When these voltages are increased the 
potential barrier between neighboring QDs is lowered and the electron wave-functions 
overlap, see two right-most QDs in Fig. 4.2. This leads (with high accuracy) to the 
effective spin-spin Heisenberg interaction JSjSj+i with J ^ 4t^/f/, where t is the 
tunneling matrix element determined by the overlap of the electron wave-functions in 
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neighboring QDs and U is the Coulomb on-site repulsion energy. The corresponding 
"clock time" can be estimated as follows: r^ ~ 1/J {h= 1). 

For the error-correcting codes to be applicable Heisenberg spin-spin interaction 
J Sj Sj+i has to be valid with the accuracy of 10"*^. Besides, to avoid correlated 
errors J should be sufficiently small in the "off" state of the two-body Hamiltonian. 
One major source of errors is the mesoscopic fluctuations. They appear in both the 
tunneling matrix element t and the relative position of the energies of the last occupied 
orbitals in neighboring QDs. In Chapter 5 we find that the constraint of having small 
errors implies keeping accurate control, at the few percent level, of several electrode 
voltages. 

The other technical problem is to control the area under the two-qubit gate pulse 
with the accuracy of 10~^. Especially, taking into account the fact that the pulse 
should be adiabatic on the time scale of tens of picoseconds so that electrons are not 
excited to the higher-lying energy levels. This issue is also discussed in details in 
Chapter 5. 

Criterion ( v) : A significant progress has been made recently in the one electron 
spin readout technique [70, 71]. Kouwenhoven and coworkers have demonstrated 
single-shot read-out of single- and two-electron spin states in a QD. An electron was 
allowed to escape from the QD or not depending on its spin state (the tunneling 
was energetically allowed only for the higher energy state). The charge of the QD 
was then measured using a quantum point contact located near the QD. With this 
technique they obtained 65% measurement visibility (the corresponding fidelity was 
82.5%). In the other slightly modified scheme they obtained even better values for 
the visibility and fidelity: 80% and 90%, respectively. Thus, high quantum efficiency 
of the measurement is reached. 

The readout time in the experiment was about 0.5 ms, whereas the energy relax- 
ation time Ti ~ 0.85 ms. Although not in principle necessary, it would be useful for 
the implementation of the error- correcting codes to reduce the readout time to about 
one hundreds of Ti time. 

In conclusion, more experiments are on the way. In particular, Kouwenhoven and 
coworkers are planning to realize swapping of the spin states in adjacent quantum 
dots. Since two-body Hamiltonian is the main building block of the spin based quan- 
tum computer these experiments will provide a valuable inside into possible directions 
of the future research. 



Chapter 5 

Spin Qubits in Multi-Electron 
Quantum Dots 

5.1 Overview 

In this chapter we study the effect of mesoscopic fluctuations on the magnitude 
of errors that can occur in exchange operations on quantum dot spin-qubits. Mid-size 
double quantum dots, with an odd number of electrons in the range of a few tens 
in each dot, are investigated through the constant interaction model using realistic 
parameters. It is found that the constraint of having short pulses and small errors 
implies keeping accurate control, at the few percent level, of several electrode voltages. 
In practice, the number of independent parameters per dot that one should tune 
depends on the configuration and ranges from one to four. 

The work in this chapter was done in collaboration with Eduardo R. Mucciolo 
and Harold U. Baranger. 

5.2 Introduction 

Since the discovery that quantum algorithms can solve certain computational 
problems much more efficiently than classical ones [56, 72], attention has been devoted 
to the physical implementation of quantum computation (QC). Among the many 
proposals, there are those based on the spin of electrons in laterally confined quantum 
dots (QD) [55], which may have great potential for scalability and integration with 
current technologies. For any successful proposal, one must be able to perform single- 
and double-qubit operations much faster than the decoherence time. In fact, all 
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logical operations required for QC can be realized if these elementary operations are 
sufficiently error free [60]. 

Single qubit operations involving a single QD will likely require precise engineer- 
ing of the underlying material or control over local magnetic fields [62]; both have 
yet to be achieved in practice. Two-qubit operations, in contrast, are already within 
experimental reach. They can be performed by sending electrical pulses to modu- 
late the potential barrier between adjacent QDs. That permits direct control over 
the effective, Heisenberg-like, exchange interaction between the qubit spins, which is 
created by the overlap between the electronic wave-functions of the QDs [55]. These 
operations are important elements in forming a basic two-qubit gate such as the 
controlled-not [73] and in the propagation of quantum information through QD ar- 
rays [55]. In fact, using three QDs instead of just one to form a logical qubit would 
allow one to perform all logical operations entirely based on the exchange interac- 
tion [69]. Thus, exchange operations will likely play a major role in the realization 
of QD qubits. A quantitative understanding of errors that occur during an exchange 
operation will help in designing optimal systems. 

The first proposal for a QD spin qubit [55] relied on having a single electron in a 
very small laterally confined QD. One advantage of such a system is that the Hilbert 
space is nominally two-dimensional. Leakage from the computational space involves 
energies of order either the charging energy or the single-particle excitation energy, 
both of which are quite large in practice (~ 1 meV ~ 10 K). Working adiabatically 

- such that the inverse of the switching time is much less than the excitation energy 

- assures minimal leakage. The large excitation energy implies that pulses of tens of 
picoseconds would be both well within the adiabatic regime and below the dephasing 
time r^ (which is typically in the nanosecond range since orbital degrees of freedom 
are involved). However, in practice, it is difficult to fabricate very small tunable 
devices [74]. Moreover, one-electron QDs may offer little possibility of gate tuning 
due to their rather featureless wave functions. 

Alternatively, a qubit could be formed by the top most "valence" electron in a QD 
with an odd number of electrons [75]. In this case, electrons filling the lower energy 
states should comprise an inert shell, leaving as the only relevant degree of freedom 
the spin orientation of the valence electron. Large QDs with 100-1000 electrons, 
while much simpler to fabricate than single electron QDs, are unsuitable because 
the excitation energy is small (~ 50/ieV ~ 0.6 K), leading to leakage or excessively 
slow exchange operations. On the other hand, mid-size QDs, with 10-40 electrons, 
are sufficiently small to have substantial excitation energies, yet both reasonable to 
fabricate and tunable through plunger electrodes. For these dots, a careful analysis 
of errors is necessary. 

Perhaps the best example of an exchange operation is the swap of the spin states 
of the two qubits. For instance, it causes up-down spins to evolve to down-up. Max- 
imum entanglement between qubits occurs when half of a swap pulse takes place - 
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a square-root-of-swap operation. Several authors have treated the problem of swap 
errors in QD systems [76, 77, 78, 79, 80, 81]. A primary concern was the occurrence 
of double occupancy (when both valence electrons move into the same QD) during 
and after the swap. However, no study so far has considered another intrinsic charac- 
teristic of electronic states in multi-electron QDs, namely, their marked dependence 
on external perturbations such as electrode voltage or magnetic field. This sensitiv- 
ity gives rise to strong sample-to-sample fluctuations arising from the phase-coherent 
orbital motion [5] . These features can make the precise control of energy levels, wave 
functions, and inter-dot couplings a difficult task. 

In this chapter we study errors and error rates that can take place during the ex- 
change operation of two spin qubits based in multi-electron QDs. We consider realistic 
situations by taking into account an extra orbital level and fluctuations in level posi- 
tions and coupling matrix elements. These lead to deviations from a pre-established 
optimal swap operation point, especially when a single-particle level falls too close to 
the valence electron level. Reasons for such fluctuations can be, for instance, (i) the 
lack of a sufficient number of tuning parameters (i.e., plunger electrodes), or (ii) the 
cross-talk between the tuning electrodes. Our results set bounds on the amount of 
acceptable detuning for mid-size QD qubits. 

This chapter is organized as follows. In Section 5.3, we introduce and justify the 
model Hamiltonian. The states involved in the exchange operation are presented in 
Section 5.4, where we also discuss the pulses and the parameters involved in the ex- 
change operations. In Section 5.5 we present the results of our numerical simulations. 
We also discuss the impact of mesoscopic effects on errors and put our analysis in the 
context of actual experiments. Finally, in Section 5.6 we draw our conclusions. 



5.3 Model System 

We begin by assuming that the double QD system can be described by the Hamil- 
tonian [12] 

H = Ha + Hb + Hab, (5.1) 

where 

a = A,B, and 

Hab = ^ (tjk Oj^bka + h.c. j . (5.3) 

j,k,(T 

Here, n^jo- = Ojo-ajo- and nB,ka = ^kJ^ka are the number operators for the single- 
particle states in the QDs (named A and B), e" denotes the single-particle energy 
levels, tjk are the tunneling amplitudes between the dots, and Ua is the charging 
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energy {a =t, i and j, k run over the single-particle states). Typically, for mid- to 
large-size QDs, the charging energy is larger than the mean level spacing. 

In the literature of Coulomb blockade phenomena in closed QDs, the Hamiltonian 
in Eq. (5.2) is known as the constant interaction model. It provides an excellent 
description of many-electron QDs, being supported by both microscopic calculations 
and experimental data [5, 12, 82, 20, 7, 83, 84]. The reasoning behind its success can 
be understood from two observations. First, mid- to large-size QDs, with more than 
ten electrons, behave very much like conventional disordered metals in the diffusive 
regime. Wavelengths are sufficiently small to resolve irregularities in the confining 
and background potentials, leading to classical chaos and the absence of shell effects 
in the energy spectrum. In this case, the single-particle states obey the statistics of 
random matrices, showing complex interference patterns and resembling a random 
superposition of plane waves. This is in contrast with the case of small, circularly 
symmetric, few-electron QDs, where shell effects are pronounced [75, 85]. 

Second, for realistic electron densities, the QD linear size is larger than the screen- 
ing length of the Coulomb interactions. In the presence of random plane waves, the 
screened interaction can then be broken up into a leading electrostatic contribu- 
tion characterized by the QD capacitance plus weak inter-particle residual interac- 
tions [12, 84]. This description becomes more accurate as the number of electrons 
gets larger since the residual interactions become weaker. The electron bunching is 
reduced as wave functions become more uniformly extended over the QD. Also, the 
increase in the number of oscillations in the wave functions leads to a self-averaging of 
the residual interactions. In this limit, one arrives at the so-called "universal Hamil- 
tonian" for QDs, containing only single-particle levels, the charging energy, and a 
mean-field exchange term [5, 12]. This Hamiltonian can be derived explicitly via a 
random-phase approximation treatment of the Coulomb interaction and the use of 
random-matrix wave functions [12, 84]. 

According to these arguments, interaction effects beyond the charging energy 
term are omitted in Eq. (5.2). In addition, the intra-dot exchange interaction, which 
tends to spin polarize the QD, is also neglected. The reason for that is the following. 
One can show that the intra-dot exchange term only affect states where there is double 
occupancy of a level. Thus, the exchange interaction constant always appears side- 
by-side with the charging energy. But in multi-electron dots, the exchange energy 
(which is at most of order of the mean level separation) is much smaller than the 
charging energy. Thus, intra-dot exchange effects are strongly suppressed by the 
charging energy. We have verified that their inclusion does not modify appreciably 
our final results. We expect the exchange interaction to become important for two- 
qubit operations only in the case of small QDs with only a few electrons, when all 
energy scales (including the mean level spacing) are of the same order. 

Thus, the simple picture where single-particle states are filled according to the 
Pauli principle up to the top most level is an appropriate description of multi-electron 
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dots [5, 12, 82, 20, 7, 83, 84]. In order to define the spin-^ qubits, both QDs should 
contain an odd number of electrons (say, 2Na—1 and 2Nb—1). The QD spin proper- 
ties are then dictated by the lone, valence electron occupying the highest level. The 
remaining electrons form an inert core, provided that operations are kept sufficiently 
slow so as not to cause particle-hole excitations to other levels. 

Experimentally, the two-qubit exchange operations also require the capability of 
isolating the QDs from each other, so that a direct product state can be prepared, 
such as 

\t) = \NAA)A(S)\NBA)B, (5.4) 

where the kets represent only the spin of the valence electron on each QD. 

5.4 Errors in Exchange Operations 

We focus our study on errors that appear after a full swap operation (which should 
result in no entanglement). Although it could in principle seem more sensible to look 
at the square-root-of-swap operation (which creates entanglement and is therefore a 
building block of logical gates), error magnitudes for the latter are straightforwardly 
related to those of the full swap operation, as we will show. We leave the discussion 
of the square root of swap to Section 5.5. 

The ideal full swap operation exchanges the valence electrons of the QD system. 
For instance, it takes the product state \i) into 

I/) = f/sw \^) = \Na, Da ® \Nb, T)b. (5.5) 

Physically, the full swap can be implemented by starting with isolated QDs, turn- 
ing on the inter-dot coupling for a time T (the pulse duration), and then turning 
it off, isolating the QDs again. For weakly coupled QDs (|t| ^ U,6e), one finds 
T ^ (7r/4) f//|tp, where t and U here represent typical values for the coupling ma- 
trix element and the charging energy, respectively (throughout we assume h = 1). 
To quantify the amount of error that takes place during the operation, we use the 
probability of not reaching |/) asymptotically, namely, 

e = l-|(/|^(+oo))r. (5.6) 

We solved numerically the time-dependent Schrodinger equation that derives from 
Eq. (5.1) for the particular but nevertheless realistic case shown in Fig. 5.1. We 
assumed that voltage tuning allows one to place the top most electron of QD A into 
an isolated single-particle state of energy e;^^ = aligned with the energy of the 
top most electron in QD B, e^ = 0. However, limited tuning ability leaves an 
adjacent empty state close in energy in QD B: efj ^^ = A. Therefore, while we 
can approximately neglect all levels but one in QD A, for QD B we needed to take 
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Figure 5.1: Schematic disposition of energy levels of a system of two QD spin qubits 
(only levels close to the top occupied state are shown). The dashed lines indicate the 
most probable transitions that can occur during the exchange operation. 



two levels into account, having hopping matrix elements denoted by ti = t^^Ng and 
^2 = tNy^,NB+i-^ To facilitate the analysis, we assumed that the dots have the same 
capacitance, C, so that Ua = Ub = U = (? jC . 

The Hamiltonian of Eq. (5.1) conserves total spin. Assuming that filled inner 
levels in both QDs are inert (forming the "vacuum" state |0)), we can span the 
5*2 = Hilbert subspace with nine two-electron basis states. According to their 
transformation properties, they can be divided into "singlet" 



lA) 
lA) 
lA) 

and "triplet" 

\Ti) 

lA) 
classes, with / = 1,2. 
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ft' -\-h' ft' 

NB+i,rNB^ ^ "NB+i,rNBi 



|0), 



(5.7) 



(5.8) 



^Though we assume time-reversal symmetry (real ti and 12), we have checked that its violation 
does not modify our results. 
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The final states that correspond to an error have either double occupancy {\Dk), 
A; = 1, . . . , 5), or an electron in the (A''^ + l)-level of QD B {\S2) and |T2)).^ In 
addition, a return to the initial state is also considered an error. It is worth noticing 
the difference between our treatment of the problem and that of Ref. [80]. In our 
case, errors come mainly from either ending in the excited single-particle state after 
the operation is over (i.e., states 15*2) and IT2)), or from "no-go" defective operations. 
In Ref. [80], errors come from having double occupancy in the final state. Double 
occupancy errors can be exponentially suppressed by adiabatically switching the pulse 
on and off on time scales larger than the inverse charging energy [55, 77, 78, 80]. 
Making pulses adiabatic on the time scale of the inverse mean level spacing for multi- 
electron quantum dots is more challenging, especially because the spacings fluctuate 
strongly both from quantum dot to quantum dot and upon variation of any external 
parameter (mesoscopic fluctuations). Therefore, multi-electron quantum dots require 
extra tunability to get around such problems. 

Very small errors, below 10^^-10"^, can, in principal, be fixed by the use of error 
correction algorithms [61]. The pulses, therefore, should be sufficiently adiabatic for 
errors to remain below this threshold. We adopted the following pulse shape: 



2 V 2r 2t J 



v{t) = -itanh — — ^ tanh — z-^ , (5.9) 

2 \ 2t 2t J 

where r is the switching time. The pulse will remain both well-defined and adiabatic 
provided that T ^ r ^ max{A~^, U~^}. Notice that this pulse is equivalent to that 
adopted in Ref. [80] up to exponential accuracy, 0{e~^''^), with T ^ r. There is no 
particular reason to believe that either performs better than the other; our choice 
was dictated by technical convenience. 



5.5 Results 

We used a standard numerical method, the so-called Richardson extrapolation [86], 
to solve the Schrodinger equation for \ilj{t)). The first step in our analysis was to find 
the optimal value of T which minimized the full swap error, as defined in Eq. (5.6), 
for a given set of parameters U, ti, and r (we used A = 1 and took t2 = ti). The 
second step was to study how this minimal error depends on r. There is actually an 
optimal interval for r, since small switching times spoil adiabaticity, while large ones 
compromise the pulse shape (when T is relatively short). Empirically, we find that 
errors related to switching times become negligible once r reaches values of about 
To = 4max{A^-'^, U^^}, provided that t <^T. In what follows, we fix r > tq. 



^While such states can lead to the full spin swap (with probability equal to 1/2), they will induce 
severe detuning errors in subsequent operations. 
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Figure 5.2: Full swap error as a function of upper level detuning in quantum dot 
B. The pulse width is optimized for A = 1, r = 6, and ti^2 = 0.2. Results for differ- 
ent charging energies are shown. Interference between different quantum mechanical 
paths in the device causes a sharp minimum. 



5.5.1 Mesoscopic Effects 



Figure 5.2 shows the full swap error as a function of A when T is fixed to its 
optimal value for A = 1. Such a situation would arise experimentally if the pulse 
is optimized for a certain configuration, but a fluctuation in level spacing occurs. 
Notice the sharp increase in error as A decreases. While increasing r reduces this 
error (by making the switching more adiabatic), very small level spacings would be 
problematic, since r can not be larger than T without sacrificing pulse shape and 
effectiveness. In order to make space for an adiabatic switching time for small A, one 
would also have to increase pulse duration. This is clear in the case of f/ = 1 (see 
Fig. 5.2): Even moderate couplings, ti,2 = 0.2, lead to larger errors, which can then 
be suppressed by decreasing ti,2 by a factor of two; however, that causes a fourfold 
increase in pulse width which may be problematic in terms of decoherence. 

The dependence of errors on fluctuations in the coupling amplitude t2 is shown 
in Fig. 5.3. Again, the pulse duration used is the optimal value obtained when 
t2 = ti = 0.2. As expected, the error grows as ^2 increases. Errors related to large 
values of t2 can also be minimized by increasing the switching time, but the same 
issues raised above appear. 

Figure 5.4 presents the error for two situations involving pulses with duration 
of about T/2, corresponding to the square-root-of-swap operation. The cases shown 
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Figure 5.3: Full swap error as a function of detuning in the coupling constant t2- 
The parameters used in the pulse width optimization are the same as in Fig. 5.2. 
Results for different charging energies are shown. 
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Figure 5.4: Comparison of error resulting from a full swap operation and two con- 
secutive square root of swap operations. The error is plotted as a function of upper 
level detuning in quantum dot B. The pulse width is optimized for A = 1, ti,2 = 0.2, 
U = 4, and r = 6. The error for a single square-root-of-swap operation is also shown. 
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Figure 5.5: Probability of liaving excessively large full swap errors (percentage) as a 
function of level spacing detuning. The solid (dashed) line corresponds to aY^I'^2 = 
(0.1) a.t U = 4 and ti = 0.2. The inset shows how the probability varies for a 
fixed width pulse (T ^ 90), but different charging energies U, when there is a large 
level-position detuning: o"a = 0.2, arj = (circles) and moderate level-position and 
coupling detunings: cr^ = 0.1, (Tr2/r2 = 0.1 (squares). 



are: (i) one, and (ii) two consecutive square-root-of-swap pulses. For comparison, the 
curve corresponding to a full swap pulse is also shown. The error for the square-root- 
of-swap operation is given by Eq. (5.6) with |/) replaced by 



I/') 



1-i 



Si) 



l + i 



Ti). 



(5.10) 



One can observe from Fig. 5.4 that error rates are nearly the same after a full swap 
operation and after two consecutive square root of swap operations. This insensitivity 
of the error to the pulse duration led us to concentrate our effort on the full swap 
operations. 

In order to establish an upper bound for QD tuning accuracy, we have performed 
simulations where both A and t2 were allowed to vary. The spacing between the levels 
in QD B was taken from a Gaussian distribution centered at A = 1, with standard 
deviation a a (e^ = was kept fixed). For the coupling amplitude, we generated 
Gaussian distributed level widths r2 = 27rt2/A, with average r2 = 27rt^/A and 
standard deviation crpj- The pulses had their widths optimized for the typical case 
where A = A = 1, f/ = 4, ti^2 = 0.2, and r = 6. For fixed values of A, o"a, 
r2, and (Jra, we generated 10,000 realizations of A and r2 and each time calculated 
the error after the application of the full swap pulse. In Fig. 5.5 we show how the 
probability of having an error larger than the 10^^ threshold depends on the energy 
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level accuracy, o"a- Two cases are considered, namely, plain and limited control of 
the inter-dot coupling constant (arj = and 0.1 r2, respectively). The data indicates 
that frequent, non-correctable errors will happen if an accuracy in A of better than 
10 percent is not achieved. 



5.5.2 Relevance for Real Quantum Dots 

To make a quantitative estimate of the impact of these results, let us consider the 
double QD setup of Jeong and coworkers [87]. In their device, each QD holds about 
40 electrons and has a lithographic diameter of 180 nm (we estimate the effective 
diameter to be around 120 nm, based on the device electron density). The charging 
energy and mean level spacing of each QD are approximately 1.8 meV and 0.4 meV, 
respectively (thus U/A ^ 4.5). If we allow for a maximal inter-dot coupling of 
^1,2 ~ 0.2 A (which yields a level broadening of about 0.25 A), we find minimal full 
swap pulse widths of about 100 ps. These values match those used in Fig. 5.5. For this 
case, switching times of 10 ps would be long enough to operate in the adiabatic regime 
and also provide an efficient and well-defined pulse shape. Thus, the combined times 
should allow for 8-10 consecutive full swap gates before running into dephasing effects 
related to orbital degrees of freedom [88] . While these numbers are yet too small for 
large-scale quantum computation, they could be sufficient for the demonstration of 
QD spin qubits. Based on Fig. 5.5, we find that accuracies in r2 of about 10% would 
make operations only limited by dephasing, and not by fiuctuat ion-induced errors. 
However, as shown in the inset, even for small QDs (typically having small f//A 
ratios), the occurrence of large errors is quite frequent when level detuning is large. 

An important issue for multi-electron QDs is their strong mesoscopic, sample- 
to-sample fluctuations in energy level position and wave-function amplitudes. Our 
results so far indicate how big an effect a given change in energy or wave-function will 
produce; now we go further and discuss how mesoscopic fluctuations more generally 
affect a collection of qubits. 

In experiments, several electrodes are placed around the QD surroundings and 
their voltages are used to adjust the lateral conflning potential, the inter-dot coupling, 
and the coupling between QDs and leads. These voltages are external parameters that 
can be used to mitigate the effects of mesoscopic fluctuations by tuning energy levels 
and wave-functions to desired values. Having that in mind, our results indicate two 
different scenarios for QD qubit implementations. 

First, if one is willing to characterize each QD pair separately and have them 
operate one by one, mesoscopic fluctuations will be irrelevant. It will be possible, 
with a single parameter per QD, say, to isolate and align energy levels reasonably 
well. Errors can be further minimized by decreasing the inter-dot coupling (thus 
increasing T) . But since QDs are not microscopically identical, each pair of QDs will 
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require a different pulse shape and duration. Multi-electron QDs are tunable enough, 
easy to couple, and much easier to fabricate than one-electron dots; therefore, multi- 
electron QDs are most appropriate for this case. 

Second, if the goal is to achieve genuine scalability, one has to operate qubits in a 
similar and uniform way, utilizing a single pulse source. In this case, T and r should be 
the same for all QD pairs. Based on our results above, one should strive to maximally 
separate the top most occupied state from all other states, occupied or empty, so as to 
reduce the possibility of leakage during operations with a fixed duration. At the same 
time, it is important to reduce inter-pair cross-talk induced by capacitive coupling 
between electrodes, as well as all inter-dot couplings except between the top most 
states of each QD. One should bear in mind that not all electrodes act independently 
- in most cases a search in a multidimensional parameter space has to be carried 
out. Thus, four tuning parameters per QD may be necessary to achieve the following 
goals: (i) find isolated, single-occupied energy level (two parameters); (ii) align this 
level with the corresponding level in an adjacent QD (one parameter); (iii) control 
the inter-dot coupling (one parameter). For parameters involved in (i) and (ii), an 
accuracy of a few percent will likely be required. Finally, control over the inter-dot 
coupling parameter, (iii), must allow for the application of smooth pulse shapes in the 
picosecond range. Our simulations also show that the pulse width must be controlled 
within at least 0.5% accuracy. Although these requirements seem quite stringent, 
recent experiments indicate that they could be met [89]. 



5.6 Conclusions 

In summary, our analysis indicate that mid-size QDs, with ten to a few tens of 
electrons, while not allowing for extremely fast gates, are still good candidates for 
spin-qubits. They offer the advantage of being simpler to fabricate and manipulate, 
but at the same time require accurate, simultaneous control of several parameters. 
Errors related to detuning and sample-to-sample fluctuations can be large, but can be 
kept a secondary concern with respect to dephasing effects provided that a sufficient 
number of independent electrodes or tuning parameters exists. 



Chapter 6 

Time Evolution of the Reduced 
Density Matrix 

6.1 Introduction 



One of the most fundamental problems in theoretical physics is that of a quantum 
mechanical system coupled to the bosonic (or fermionic) bath [90, 91]. One such 
example is presented in the next chapter where we consider double QD charge qubit 
coupled to the acoustic phonon bath. Formal solution to the problem can be found 
[see Eqs. (6.17) and (6.18) below], nonetheless to do any practical calculation one has 
to (a) make simplifying assumptions like Born and Markov approximations in the 
derivation of the Redfield equation [92, 93] or (b) implement sophisticated numerical 
integration procedure like QUAPl (quasiadiabatic propagator path integral) [94, 95, 
96, 97, 98]. 

Time evolution of the isolated quantum mechanical system is governed by the 
Schrodinger equation. Therefore, (i) its energy is conserved and (ii) its phase coher- 
ence is preserved over time. Coupling to the bath leads to the energy relaxation - 
decoherence and loss of phase coherence - dephasing. These processes are important 
for quantum computations because they may potentially destroy the phase coher- 
ence of both unentangled single qubit states and especially entangled states of many 
qubits. 

In our treatment of the problem we choose to go along the lines of Born and 
Markov approximations which eventually lead us to the Redfield equation [92, 93] 
because (i) we are interested in time evolution of the reduced density matrix at 
sufficiently long times and (ii) this formalism allows us to obtain analytical expressions 
for both qubit decoherence and dephasing rates. These analytical expressions contain 
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valuable information about the dependence of the rates on different parameters of the 
double QD setup. 

In this chapter we give formal derivation of the Redfield equation making simpli- 
fying assumptions along the way. We employ the projection operator technique [99] 
as it provides the most transparent way of the derivation. 

The plan of this chapter is as follows. In Section 6.2 we formulate the problem. 
In Section 6.3 we find formal solution of the problem. In Section 6.4 we assume that 
the quantum mechanical system and the bath are unentangled at t = and also 
make Born approximation. Finally, in Section 6.5 we make Markov approximation 
and finalize the derivation of the Redfield equation. 

6.2 Formulation of the Problem 

The Hamiltonian of the total system is 

H{t) = Hs{t) + HB + V, (6.1) 

where Hs{t) is the Hamiltonian of the system, Hb is the Hamiltonian of the bath, 
and V is the Hamiltonian of their interaction. We assume that it can be factorized 
as follows: 

V = J2K^<^^, (6-2) 

i 

where {Ki} are the operators in the system's of interest Hilbert space and {^i} are 
the operators in the bath Hilbert space. Time evolution of the total system obeys 
Schrodinger equation for the density matrix: 

tj^p{t) = [H{t),p{t)]. (6.3) 

Our goal is to find similar equation for the time evolution of the reduced (system's of 
interest) density matrix o"(t). 

6.3 General Solution 

Let us write p(t) as follows: 

pit) = fiHB)ait)+vit), (6.4) 

where 
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One can always do this because no assumption is made about i]{t). Now let us 
substitute Eq. (6.4) into Eq. (6.3). We obtain 

t \f{HB)a{t)+f,] = fiHB)[Hsit),ait)] + [Hsit),r]it)] (6.6) 

+ [HB,v{t)] + [VJXHBMt)] + [V,vit)]. 

To find the equation on a{t) let us trace left and right hand sides of Eq. (6.6) over 
bath degrees of freedom. The following identities are useful in simplifying the traces: 

tnr]{t) = 0, tnf{HB) = l, tn[Vf{HB)] = J2^it^b[^^f(HB)] = 0. (6.7) 

i 

The last identity should be valid for the most practical applications (it is certainly 
valid for the phonons coupled to the electric charge, for example). Then the equation 
for the reduced density matrix becomes 

t&{t) = [Hs{t),a{t)]+tn[V,r]{t)], (6.8) 

where 

tr,[V,i]it)] = J2[K^,tn{<!>,l^(t)}]. (6.9) 

i 

However there is one problem: Eq. (6.8) contains ri{t) as well. Therefore, one has 
to find the second (coupled to the first one) equation on i]{t). This can be done by 
acting with the projection operator P = 1 — f{HB)trb on both sides of Eq. (6.6). A 
little bit more involved calculation leads to the following result: 

tf,{t) = [iJ5(t),r7(t)] + [i/B,r7(t)] + [\/,77(t)] + [\/,/(iJB)^(t)]-/(iJB)tr4\/,r7(t)]. (6.10) 

One can write Eqs. (6.8) and (6.10) in a more compact form as follows: 

a{t) = -i'Hs{t)(y{t) -itnVr]{t) (6.11) 

m = -^ [Hsit) +nB + PVUt) - t VfiHBMt), (6.12) 

where Hs, 'Hb, and V are the Liouvillian operators: 'Hs(t)a{t) = [Hs{t),a{t)], 
Vv{t) = [V,r]it)],etc. 

Our further plan is as follows. Firstly, let us solve Eq. (6.12), that is, find the 
solution for ri{t) - it is going to have some functional dependence on a{t). Secondly, 
let us substitute the solution in Eq. (6.11) and simplify the expression obtained. 

Eqs. (6.11) and (6.12) have the following form: 

r]{t) = A{t)r]{t) + B{t) (6.13) 

r]{0) = r]o. (6.14) 
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This differential equation has a standard solution which can be written as follows 

r]{t) = S{t,0)r]o+ dt'S{t,t')B{t'), where 5(t, t') = te-^*''^''^^^) (6.15) 

Jo 

and T is a time-ordering operator. How do we understand S{t^ ^)flo term in Eq. (6.15), 
for example? It should be understood as follows: 

5(t,0)r7o = fe-^o'^"'4(")r7o = f e-^o ^"^(") ^gf^e' J'o '^-^W, (6.16) 

where T_ is an anti-time-ordering operator. Eq. (6.15) solves Eq. (6.12). Let us 
substitute the solution in Eq. (6.11). Thus, we obtain closed form differential equation 
for a{t): 

a{t) = -i ns{t)a{t) - i tnVSit, 0)rio - tnV f dt'S{t, t')Vf{HB)(r{t'), (6.17) 

Jo 

where 

5(t,t') =Te-^-^*''^"[^^(")+^^+^^l. (6.18) 

An important note is due at this time. At t = 0: p(0) = f{H 3)1^0 + ^o- Therefore, if 
system and bath are not entangled at t = then r^o = and second term in Eq. (6.17) 
simply vanishes. 



6.4 770 = and Born Approximations 

Eq. (6.17) may look simple but it is difficult to solve mainly due to complex 
structure of the S matrix. Let us make a couple of simplifying assumptions. Firstly, 
let us assume that r^o = then second term in Eq. (6.17) vanishes. Secondly, let us 
assume that the system is weakly coupled to the bath {V is small) and keep only the 
leading term in powers of V expansion of the third term (Born approximation): 

S{t,t') = 5o(t,t') + 0(V), where 5o(t,t') = f e-'-^*''^"[^«(")+^«l. (6.19) 

Thus, in r^o = and Born approximations, Eq. (6.17) becomes 

&{t) = -ins{t)a{t) - tnV f dt'Soit,t')VfiHBMt'), (6.20) 

Jo 

where the first term on the right hand side is the conventional Liouvillian term cor- 
responding to the time reversible dynamics and the second term is a dissipative term 
(or time irreversible part). 

The dissipative term can be further simplified. First, let us simplify Sq: 

So{t,t') = Ss{t,t')e-'^'-''^^^ , where Ss{t,t') = fe-'^^'"'''^'^'=\ (6.21) 



Then, the dissipative term becomes 

- / dt'tn \v,Ss{t,t') [e-*(*-*')^^\/e'(*-*')^^,/(i/B)a(tOll • (6.22) 

Further calculations should be performed in the following order: (i) substitute V = 
Yli^i^i'j (ii) i^^^ advantage of tvb - define bath correlation functions: 

B,,{t)=tn[<^,{t)^jf{H,)], (6.23) 

B,,i-t)=tn[<!>,<!>,{t)f{H,)], (6.24) 

[B,,i-t)]^ = B,,ity, (6.25) 

(iii) to transform a{t — y) to a{t) use:^ 

Ssit, t - y) ait -y) = a{t) + O(V'). (6.26) 

The final form of the equation for the time evolution of the reduced density matrix 
in ?7o = and Born approximations is 

&{t) = -t [Hs{t), ait)] + Yl t^^^^) ^W' ^^l + H-c, (6.27) 

j 

where 

^j(t) = y^ [ dyB^,{y)Ss{t,t-y)K, (6.28) 

^ Jo 

= V / rfy%(|/)fe''-^*%''"^^(")i^,f_e*-^'^''"^^("\ (6.29) 

i Jo 

The structure of A-matrices is quite complex. Moreover, at each time slice in the 
integration of Eq. (6.27) one has to numerically calculate A-matrices all over again. 

6.5 Markov Approximation on Top: Redfield Equa- 
tion 

One conventional way to simplify Eq. (6.29) and make Eq. (6.27) local in time 
is to assume that bath has no memory. That is, bath correlation time Tc is the 
smallest time scale in the problem. This is called Markov approximation. In this 
approximation one can only find evolution of the system at sufficiently long times: 



^Eq. (6.26) implies the validity of the Markovian approximation as well. 
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t > Tc- Since Bji{y) -^ for y > Tc, only small values of y contribute in Eq. (6.29). 
Hence, 

/ dx Hs{x) ^ y Hs{t) (6.30) 

Jt-y 

and upper limit of the integration over y in Eq. (6.29) can be extended to oo. Thus, 
A-matrices in Born-Markov approximation can be written as follows: 

Aj(t) = V/ rfyi?,,(i/)e-^^^«Wi^,e^^^^W. (6.31) 

i Jo 

Eq. (6.27) with A-matrices given by Eq. (6.31) is local in time and called Redfield 
equation [92, 93, 99, 100]. 



Chapter 7 

Phonon Decoherence of a Double 
Quantum Dot Charge Qubit 

7.1 Overview 

In this chapter we study decoherence of a quantum dot charge qubit due to 
couphng to piezoelectric acoustic phonons in the Born-Markov approximation. After 
including appropriate form factors, we find that phonon decoherence rates are one 
to two orders of magnitude weaker than was previously predicted. We calculate the 
dependence of the Q-factor on lattice temperature, quantum dot size, and interdot 
coupling. Our results suggest that mechanisms other than phonon decoherence play 
a more significant role in current experimental setups. 

The work in this chapter was done in collaboration with Eduardo R. Mucciolo 
and Harold U. Bar anger. 

7.2 Introduction 

Since the discovery that quantum algorithms can solve certain computational 
problems much more efficiently than classical ones [101], attention has been devoted 
to the physical implementation of quantum computation. Among the many proposals, 
there are those based on the electron spin [55, 69] or charge [102, 103, 104, 105, 81, 106] 
in laterally confined quantum dots, which may have great potential for scalability and 
integration within current technologies. 

Single qubit operations involving the spin of an electron in a quantum dot will 
likely require precise engineering of the underlying material or control over local 
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Figure 7.1: Schematic Coulomb blockade stability diagram for a double quantum dot 
system at zero bias [107]. (A^i, N2) denotes the number of excess electrons in the dots 
for given values of the gate voltages Vi and V2. The solid lines indicate transitions in 
the total charge, while the dotted lines indicate transitions where charge only moves 
between dots. The point A marks the qubit working point. 



magnetic fields [62]; both have yet to be achieved in practice. In contrast, single qubit 
operations involving charge in a double quantum dot (DQD) [107] are already within 
experimental reach [88, 108, 109]. They can be performed either by sending electrical 
pulses to modulate the potential barrier between the dots (tunnel pulsing) [104, 105, 
106] or by changing the relative position of the energy levels (bias pulsing) [88, 108]. 
In both cases one acts on the overlap between the electronic wave functions of the 
dots. This permits direct control over the two low-energy charge states of the system - 
the basis states |1) and |2) of a qubit: Calling A^i (A^2) the number of excess electrons 
in the left (right) dot, we have that |1) = (1,0) and |2) = (0, 1). 

The proposed DQD charge qubit relies on having two lateral quantum dots tuned 
to the (1,0) ^^ (0, 1) transition line of the Coulomb blockade stability diagram (see 
Fig. 7.1). Along this line, an electron can move between the dots with no charging 
energy cost. An advantage of this system is that the Hilbert space is two-dimensional, 
even at moderate temperatures, since single-particle excitations do not alter the 
charge configuration. Leakage from the computational space involves energies of 
order the charging energy which is quite large in practice (~ 1 meV ~ 10 K). In the 
case of tunnel pulsing, working adiabatically - such that the inverse of the switching 
time is much less than the charging energy - assures minimal leakage. The large 
charging energy implies that pulses as short as tens to hundreds of picoseconds would 
be well within the adiabatic regime. However, the drawback of using charge to build 
qubits is the high decoherence rates when compared to spin. Since for any successful 
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qubit one must be able to perform single- and double-qubit operations much faster 
than the decoherence time, a quantitative understanding of decoherence mechanisms 
in a DQD is essential. 

In this chapter, we carry out an analysis of phonon decoherence in a DQD charge 
qubit. During qubit operations, the electron charge movement induces phonon cre- 
ation and annihilation, thus leading to energy relaxation and decoherence. In order to 
quantify these effects, we follow the time dependence of the system's reduced density 
matrix, after tracing out the phonon bath, using the Redfield formalism in the Born 
and Markov approximations [99, 100]. 

Our results show that decoherence rates for this situation are one to two orders 
of magnitude weaker than previously estimated. The discrepancy arises mainly due 
to the use of different spectral functions. Our model incorporates realistic geomet- 
ric features which were lacking in previous calculations. When compared to recent 
experimental results, our calculations indicate that phonons are likely not the main 
source of decoherence in current DQD setups. 

The chapter is organized as follows. In Section 7.3, we introduce the model used 
to describe the DQD, discuss the coupling to phonons, and establish the Markov 
formulation used to solve for the reduced density matrix. In Section 7.4 we study de- 
coherence in a single-qubit operation, while in Section 7.5 we simulate the bias pulsing 
experiment of Refs. [88, 108]. Finally, in Section 7.6 we present our conclusions. 



7.3 Model System 

We begin by assuming that the DQD is isolated from the leads. The DQD and 
the phonon bath combined can then be described by the total Hamiltonian [81] 

H = Hs + Hb + Hsb, (7.1) 

where Hs and Hb are individual DQD and phonon Hamiltonians, respectively, and 
HsB is the electron-phonon interaction. We assume that gate voltages are tuned 
to bring the system near the degeneracy point A (Fig. 7.1) where a single electron 
may move between the two dots with little charging energy cost. To simplify the 
presentation, only one quantum level on each dot is included; -Ei(2) denotes the en- 
ergy of an excess electron on the left (right) QD (possibly including some charging 
energy). Likewise, spin effects are neglected.^ Thus, in the basis {|1), |2)}, the DQD 



^In the situation we envision, the total number of electrons would be odd: The ground state of 
each dot in the absence of the excess electron would have 5* = 0, so that the dot with the excess 
electron would have spin half. Other situations are, of course, possible, such as the total number of 
electrons being even with S — for both dots when the qubit is in state |1) and S = 1/2 for both 
in state |2); in the latter case, there would be a singlet and triplet state of the DQD with a small 
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Figure 7.2: Geometry of the double quantum dot charge qubit. 



Hamiltonian reads 






(7.2) 



where az^x are Pauh matrices, e{t) =Ei—E2 is the energy level difference, and v{t) is 
the tunneling amplitude connecting the dots. Notice that both e and v may be time 
dependent. The phonon bath Hamiltonian has the usual form {h = 1) 



HB = ^UJq 6]j6q, 



(7.3) 



where the dispersion relation tUq is specified below. The electron-phonon interaction 
has the linear coupling form [81, 110], 



^5^ = EE"?^^ W+M 



(7.4) 



q i=l 



where Nj is the number of excess electrons in the i-th dot and a 



(i) 



X g-»qR. p. 



tq — ^vq"^ - J- ^{^.)1 

with Ri = and R2 = d the dot position vectors, see Fig. 7.2. The dependence of the 
coupling constant Aq on the material parameters and on the wave vector q will be 
specified below. The dot form factor is 



^»(q) 



rfVn,(r)e-^^-' 



(7.5) 



exchange splitting. Such comphcations do not effect the underlying physics that we discuss, and so 
we neglect them. 
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where nj(r) is the excess charge density in the i-th dot. With no significant loss of 
generahty, we will assume that the form factor is identical for both dots and, therefore, 
drop the i index hereafter. In the basis {|1), |2)}, after dropping irrelevant constant 
terms, the electron-phonon interaction simplifies to 

HsB = K<l>, (7.6) 

where 

K=-(T, and $ = 5^ ^q (&i + fe-q) , (7.7) 

q 
with (yfq = AqP(q) (l — e~*'^*^). The phonons propagate in three dimensions, while 
the electrons are confined to the plane of the underlying two-dimensional electron 
gas (2DEG). Notice that the electron-phonon coupling is not isotropic for the DQD 
(Fig. 7.2): Phonons propagating along = and any 9 do not cause any relaxation, 
while coupling is maximal along (j) = 6 = 7r/2 direction. We neglect any mismatch in 
phonon velocities at the GaAs/AlGaAs interface, where the 2DEG is located. 

We now proceed with the Born-Markov-Redfield treatment of this system, see 
Chapter 6. While the Born approximation is clearly justified for weak electron-phonon 
interaction, the Markov approximation requires, in addition, that the bath correlation 
time is the smallest time scale in the problem. These conditions are reasonably 
satisfied for lateral GaAs quantum dots, as we will argue below. 

Let us assume that the system and the phonon bath are disentangled at t = 0. 
Using Eqs. (7.2), (7.3), and (7.6), we can write the Redfield equation for the reduced 
density matrix p{t) of the DQD (see Chapter 6), 

Pit) = -nHsit),pit)] + {[Ait)pit),K]+R.c.} . (7.8) 

The first term on the right-hand side yields the Liouvillian evolution and the other 
terms yield the relaxation caused by the phonon bath. The auxiliary matrix A is 
defined as 

A(t)=/ rfr5(r)e-^"^^Wire*"^sW (7.9) 

Jo 
where _B(r) =Trb{$(r)$(0)/(/JB)} is the bath correlation function, <l>(r) = e'^^^ $ e~''^^'^, 
and JXHb) = e^^^^ /TTi){e^^^^}, with f3 = 1/T the inverse lattice temperature 
{kB = 1). 

Using Eq. (7.3) in the definition of the bath correlation function, we find that the 
latter can be expressed in the form 



B{t)= dLuiy{iu){e'^^nBH + e-'^'^[l + nBiuj)]}, (7.10) 

Jo 

where nB{oo) is the Bose- Einstein distribution function and 
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is the spectral density of the phonon bath. 

We now speciahze to hnear, isotropic acoustic phonons: C(jq = s|q|, where s is the 
phonon velocity. Moreover, we only consider coupling to longitudinal piezoelectric 
phonons, neglecting the deformation potential contribution. For bulk GaAs, this is 
justifiable at temperatures below approximately 10 K (see Appendix B). Thus, 



IX |2 



n\q\ '' 



(7.12) 



where (7ph is the piezoelectric constant in dimensionless form ((7ph~0.05 for GaAs [110, 
111]) and Q is the unit cell volume. 

The excess charge distribution in the dots is assumed Gaussian: 



n[r} 



5{z) 



27ra2 



exp 



2a2 



(7.13) 



This is certainly a good approximation for small dots with few electrons, but becomes 
less accurate for large dots. The resulting form factor reads 



P(q) 



-{g2+g2)aV2 



(7.14) 



Note that this expression differs from that in Refs. [104, 105, 106] where a three- 
dimensional Gaussian charge density was assumed. 

Using Eqs. (7.12) and (7.14), as well as the DQD geometry of Fig. 7.2, we get 



uiu 



fi'ph^ 



tt/2 



dO sin 9 exp 







uj^a" 



sm^O 



l-Ji 



ood 



sin 9 



(7.15) 



It is instructive to inspect the asymptotic limits of this equation. At low frequencies, 
v{uj -^ 0) ~(7ph (i^cj^/Gs^; thus, the phonon bath is superohmic. At high frequencies, 

o2 



iy{uj ^ oo) ^ -^ — / - 



where 



a 



dx X e 



Jn 



d 

—X 

a 



(7.16) 



(7.17) 



Notice that the spectral function does not have the exponential decay familiar from 
the spin-boson model, but rather falls off much more slowly: z/(c<j^ oo) occj^^. This 
should be contrasted with the phenomeno logical expressions used in Ref. [81]. 

The characteristic frequency of the maximum in /^(o;) is r~^ = s/a. For typical 
experimental setups, a~50 nm while s^5 x 10^ m/s for GaAs, yielding Tc ~ 10 ps 
{t^^ ^65/ieV). Thus, the Markovian approximation can be justified for time scales 
t > Tc and if all pulse operations are kept adiabatic on the scale of Tc. 
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7.4 Decay of Charge Oscillations 

One can operate this charge qubit in two different ways: (i) by pulsing the tunnel- 
ing amplitude v{t) keeping e constant, or (ii) by changing the energy level difference 
e{t) keeping v constant (bias pulsing). Tunnel pulsing seems advantageous as it im- 
plies fewer decoherence channels and less leakage. However, a recent experiment used 
a bias pulsing scheme [88, 108]. 

Our system's Hilbert space is two-dimensional by construction [see Eq. (7.2)], 
hence there is no leakage to states outside the computational basis. We can, therefore, 
use square pulses instead of smooth, adiabatic ones. This not only allows us to 
analytically solve for the time evolution of the reduced density matrix, Eq. (7.8), but 
also renders our results applicable to both tunnel and bias pulsing. Indeed, in both 
regimes one has e(t) =0 and v(t) = t>m for t > 0, taking that the pulse starts at 
t = 0. Let us assume that the excess electron is initially in the left dot: pii(O) = 1 
and pi2(0) = 0. In this case, since the coefficients on the right-hand side of (7.8) are 
all constants at t > 0, we can solve the Redfield equation exactly (see Appendix C 
for details). As p{t) has only three real independent components, the solution is 



Repuit) 
Impi2(t) 



1 1 71 Tl 

— I — e 2"*(coscjt H sincjt), 

2 2 ^ 2u '' 



-^(l-e--^')tanh^, 

2Wm + 72 -TLi . 

■ e 2 sm a;t, 



where 



71 

72 



2uj 






1/2 



n 



Vr, 



2 v{2vm)coi\i-^ , 



y 



^— YZ/(2t;„,t/)coth^. 



(7.18) 
(7.19) 
(7.20) 

(7.21) 
(7.22) 
(7.23) 



Note that 7i,2<^'i^m- We extract the customary energy and phase relaxation times, 
Tl and T2, by rotating to the energy eigenbasis {|— ), |+)}: 



P—{t) = --Repi2(t), 

P-+{t) = --+Pn{t)+ilmpi2{t). 



(7.24) 
(7.25) 



Then, the damping of the oscillations in the diagonal matrix elements is the signature 
of energy relaxation, while the phonon-induced decoherence is seen in the exponential 
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decay of the off-diagonal elements. For the DQD, we find Ti =7^^ 
the decoherence time. 



The quality factor of the charge oscillations in Eq. (7.1^ 
Eqs. (7.21), (7.22), and (7.15), we find that 



is Q 



Q 



4tanh(t>m/T) 



TT^^ph 



dx 



\^l — X 



'(Vm/^a)'^ X 



Jn 



dVr, 



a Un 



and T2 = 27^^ ^ for 
uj/n'-^i. Using 

-1 

, (7.26) 



Vm, lattice 



where uja = s/2a. The Q-factor depends on the tunneling amplitude 
temperature T, dot radius a, and interdot distance d. 

Several experimental realizations of DQD systems recently appeared in the liter- 
ature [88, 108, 109, 87, 89]. In principle, all these setups could be driven by tunnel 
pulsing to manipulate charge and perform single-qubit operations. To understand 
how the Q-factor depends on the tunneling amplitude Vm in realistic conditions, let 
us consider the DQD setup of Jeong and coworkers [87]. In their device, each dot holds 
about 40 electrons and has a lithographic diameter of 180 nm. The effective radius 
a is estimated to be around 60 nm based on the device electron density. Therefore, 
d/a^3. The lattice base temperature is 15 mK. Introducing these parameters into 
Eq. (7.26), one can plot Q-factor as a function of Vm or, equivalently, as a function of 
the period of the charge oscillations P = 2TT/ujK.TT/vm- This is shown in Fig. 7.3. 

To stay in the tunnel regime Vm should be smaller than the mean level spacing 
of each QD, approximately 400 /leV in the experiment [87]. Therefore, in Fig. 7.3 
we only show the curve for Vm up to 100 /leV. One has to recall that at these values 
the Markov approximation used in the Redfield formulation is not accurate (see end 
of Section 7.3), and so our results are only an estimate for Q. For strong tunneling 
amplitudes, when 25/ieV < Vm < lOOyueV, the largest value we find for Q is close to 
100. For weak tunneling with Vm < 25 /xeV, the situation is more favorable and larger 
quality factors (thus relatively less decoherence) can be achieved. Nevertheless, the 
one-qubit operation time, which is proportional to the period, grows linearly with Q 
in the region of Vm —^ 0, as shown in the inset of Fig. 7.3. Therefore, at a certain 
point other decoherence mechanisms are going to impose an upper bound on Q. 

The minimum of Q in Fig. 7.3 occurs when Vm coincides with the frequency at 
which the phonon spectral density is maximum. It corresponds to the energy splitting 
between bonding and anti-bonding states of the DQD, 2vm, being approximately equal 
to the frequency of the strongest phonon mode s/a: Vm — uJa- 

From Fig. 7.3, it is evident that one can reach certain values for the Q-factor (say, 
Q = 100) at both weak {vm ^ 4.6 /leV ~ 53 mK) and strong (f^ ~ 93 /xeV ~ 1.1 K) 
tunneling. However, these two regimes are not equally convenient. From Eq. (7.26), 
it is clear that the temperature dependence of the Q-factor is fully determined by the 
bonding-antibonding splitting energy 2vm'- Q(T) =Q(0) tanh(fm/T). We notice that 
Q(T) ^ (5(0) if T <^ Vm] therefore, the Q-factor is less susceptible to temperature 
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Figure 7.3: The charge oscillation Q-factor as a function of the tunneling amplitude 
Vm (lower scale) and of the oscillation period P (upper scale) for a GaAs double 
quantum dot system. The lattice temperature is 15 mK and the dot radius and 
interdot distance are 60 nm and 180 nm, respectively. The inset shows the relation 
between P and Q at small tunneling amplitudes (large periods). 
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Figure 7.4: The charge oscillation Q-factor as a function of the lattice temperature. 
Inset: as a function of the dot radius for a fixed ratio dja = 3. The solid (dashed) line 
corresponds to the weak Vm — 53 mK (strong fm, — 1-1 K) tunneling regime. Other 
parameter values are equal to those in Fig. 7.3. 
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variations for strong tunneling (Fig. 7.4). Another parameter that influences the Q- 
factor is the dot radius, which controls the frequency of the strongest phonon mode, 
s/a. In the strong tunneling regime (dashed curve in the inset to Fig. 7.4), one has 
to increase the QD size to improve the Q-factor. This would reduce the energy level 
spacing, hence only moderate improvement in Q- factor is possible. In contrast, in 
the weak tunneling regime (solid curve in the inset to Fig. 7.4) one has to reduce 
the QD size. This can lead to a significant (up to one order of magnitude) Q-factor 
improvement. 



7.5 Bias Pulsing 

In a recent experiment [88, 108], Hayashi and coworkers studied charge oscillations 
in a bias-pulsed DQD. In this regime the energy difference between the left and right- 
dot single-particle energy levels is a function of time: e{t) = eQu{t). A typical profile 
used for pulsing is 

1 /^ , t + W/2 , t~W/2\ ,^ ^, 

u{t) = 1 - - ( tanh ^^' - tanh ^^' \ , (7.27) 

where W represents the pulse width and r controls the rise and drop times. During 
bias pulsing, the tunneling amplitude is kept constant. In Refs. [88, 108], the difference 
in energy levels was induced by applying a bias voltage between left and right leads 
(and not by gating the dots separately). For their setup, the maximum level splitting 
amplitude was Eq ~ 30/ieV and r ~ 15 ps, corresponding to an effective ramping 
time of about 100 ps.^ The tunneling amplitude was kept constant and estimated as 
f ~ 5 /ieV, which amounts to charge oscillations with period P ?^ 1 ns. The lattice 
temperature was 20 mK. Each quantum dot contained about 25 electrons and the 
effective dot radius is estimated to be around 50 nm based on the device electron 
density. From the electron micrograph of the device one finds d ~ 225 nm, hence 
ci/a~ 4.5. When substituting these values into Eq. (7.26), one finds (5~54. 

However, from the experimental data one observes Q ~ 3. Low Q-factors were also 
obtained by Petta and coworkers in an experiment where coherent charge oscillations 
in a DQD were detected upon exciting the system with microwave radiation [109]. 
Other mechanisms of decoherence do exist in these systems, such as background 
charge fluctuations [113] and electromagnetic noise emerging from the gate voltages. 
Our results combined with the recent experiments indicate that these other mecha- 
nisms are more relevant than phonons. 

We now turn to yet another possible source of decoherence: Leakage to the leads 



^The value of tw 15 ps is obtained by fitting Eq. (7.27) to the experimental pulse with an effective 
ramping time of 100 ps, Ref. [112]. 
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when the pulse is on.'^ To illustrate this alternative source of damping of charge os- 
cillations, we simulate the bias-pulsing experiment of Refs. [88, 108] by implementing 
a rate equation formalism similar to that used in Ref. [112]. The formalism is based 
on a transport theory put forward for the strongly biased limit [115, 116]. First, we 
find the stationary current Iq through the DQD structure when the pulse is off (that 
is, the bias is applied) [116]: 

' ''"" "' ^. (-.28) 



L + ^R -y2 _^ El£r _^ 



(rL+rfl)2 



where e is the elementary charge. ^l{r) is the partial width of the energy level in 
the left (right) dot due to coupling to the left (right) lead (when the bias is applied); 
in the experiment [88, 108], T^r ^ 30/ieV. On the other hand, when the pulse is 
on, the stationary current is zero. We now apply the pulse e{t) and measure the 
current I{t). In the experiments, the level widths T^fi decrease upon biasing the 
system. To include that effect here, we also pulse them: Tii^t) =7l + (F^ — jL)u(t) 
and analogously for F/j, where •jLiR) is the residual leakage to the left (right) lead 
when the pulse is on. We use 'yL,R = 0.3 /xeV, even though the real leakage in the 
experiment was likely much smaller. To obtain the response current one subtracts 
the stationary component: Ircsp(t)=I{t) — Iou(t). 

Figure 7.5(a) shows the response current for a pulse of width W^ = 4 ns and 
r = 30 ps. The latter is approximately twice as large as in the experiment and is 
chosen to enhance the effect. In Refs. [88, 108], pulses were applied at a frequency 
/ = 100 MHz. The average number of electrons transfered from the left to the right 
lead per cycle minus that in the stationary regime is [112] 

/■I// 
n= / rft /resp(t)/e. (7.29) 



(In the simulations there is no need to apply a sequence of pulses.) Notice that n 
oscillates as a function of the pulse width W [see Fig. 7.5(b)] as observed in the ex- 
periment. Two main conclusions can be drawn from our simulation. First, the larger 
r, the smaller the visibility of the charge oscillations [112]. Second, the larger the 
leakage rates '~fL,R when the pulse is on, the stronger the damping of the oscillations. 
While the damping due to leakage is presumably too weak an effect to discern in the 
data presented in Refs. [88, 108], the loss of visibility due to finite r is likely one of 
the causes of the small amplitude seen experimentally. 



^For another source of dccoherence due to coupling to the leads, see Ref. [114]. 
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Figure 7.5: (a) The response current Jresp(t)/e in ns^^ as a function of time for a 
pulse with W^ = 4 ns and r = 30 ps. (b) Number of electrons transferee! between left 
and right leads, as defined in Eq. (7.29), as a function of the pulse width W. 
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7.6 Conclusions 

The main conclusion of the chapter is that, under reahstic conditions, phonon de- 
coherence is one to two orders of magnitude weaker than expected [81, 104, 105, 106]. 
The analytical expression for the Q-factor given in Eq. (7.26) was found using an ex- 
pression for the phonon spectral density, Eq. (7.15), which takes into account impor- 
tant information concerning the geometry of the double quantum dot system. In a pre- 
vious work [81] an approximate, phenomenological expression, /^(u;) oc u;exp(— cj/cJc), 
was utilized in the treatment of charge qubits. There is a striking difference be- 
tween these two expressions in both the high- and low-frequency limits. Moreover, 
an arbitrary coupling constant was adopted in Ref. [81] to model the electron-phonon 
interaction while our treatment uses a value known to describe the most relevant 
phonon coupling in GaAs. On the other hand, other previous work [104, 105, 106] 
assumed a spherically symmetric excess charge distribution in the dot while we have 
assumed a two-dimensional pancake form. These differences account for most of the 
discrepancy between the present and previous results. 

Based on these findings we conclude that phonon decoherence is too weak to 
explain the damping of the charge oscillations seen in recent experiments [88, 108, 
109]. Charge leakage to the leads during bias pulsing is an additional source of 
damping, as shown in Fig. 7.5(b); however, for realistic parameters [88, 108, 112], it 
turns out to be a weak effect as well. Hence, other decoherence mechanisms, such 
as background charge fluctuations or noise in the gate voltages, play the dominant 
role [113]. 

There are two distinct ways to operate a double quantum dot charge qubit: (i) 
by tunnel pulsing or (ii) by bias pulsing. Tunnel pulsing seems advantageous due to 
the smaller number of possible decoherence channels. In addition, the bias pulsing 
scheme, in contrast to tunnel pulsing, introduces significant loss of visibility in the 
charge oscillations. 

In this chapter we did not attempt to study leakage or loss of fidelity due to 
non-adiabatic pulsing, which are both important issues for s^m-based quantum dot 
qubits (see Chapter 5). Moreover, we have not attempted to go beyond the Markov 
approximation when deriving an equation of motion for the reduced density matrix. 
Both of these restrictions in our treatment impose some limitations on the accuracy 
of our results, especially for large tunneling amplitudes. 

Finally, it is worth mentioning that some extra insight would be gained by mea- 
suring the Q-factor as a function of the tunneling amplitude Vm experimentally. Such 
a measurement would allow one to map the spectral density of the boson modes re- 
sponsible for the decoherence. This would provide very valuable information about 
the leading decoherence mechanisms in double quantum dot systems. 



Chapter 8 



Conclusion 



8.1 General Overview 

The continuous minituarization of the integrated circuits is going to affect the 
underlying physics of the future computers. This new physics first came into play as 
the effect of Coulomb blockade in the electron transport through the small conducting 
island. Then, as the size of the island L continued to shrink further, quantum phase 
coherence length L^ became larger than L leading to the mesoscopic fluctuations 
- fluctuations of the island's quantum mechanical properties upon small external 
perturbations. Quantum coherence of the mesoscopic systems is essential for building 
reliable quantum computer. Unfortunately, one can not completely isolate the system 
from the environment and its coupling to the environment inevitably leads to the loss 
of coherence or decoherence. All these effects are to be thoroughly investigated as 
the potential of the future applications is enormous. 



8.2 Main Conclusions of My Work 

In Chapter 2 I studied Coulomb blockade oscillations of the linear conductance 
through a quantum dot weakly coupled to the leads via multichannel tunnel junctions. 
To obtain analytic results I have assumed that the energy levels in the QD are equally 
spaced. The electron-electron interactions in the QD have been described by the 
constant interaction model, though thermal excitations with all possible spins were 
taken into account. 

Firstly, I found the expression for the linear conductance in the spinless case 
[Eq. (2.22)]. It is valid at arbitrary values of Ec-, SE and T. Then, in Section 2.5, I 
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applied the spinless case theory resuh to the problem of the transport via quantum 
dot in the quantum Hall regime. The quantum dot energy levels in this case are 
equidistant with the spacing given by Eq. (2.35). 

Linear conductance in the case of spin-| electrons at arbitrary values of Ec, SE, 
and T is given by Eq. (2.55). It is plotted in Fig. 2.9 for the charging energy of the 
quantum dot equal to mean level spacing. Eq. (2.55) can be used as an experimental 
data fit to extract the charging energy and mean level spacing at given temperature. 

Though I do not expect my quantitative results to precisely describe a quantum 
dot with random energy levels, they certainly give correct order of magnitude for the 
conductance oscillations and their generic features. 

In Chapter 3 I studied corrections to the spacings between Coulomb blockade 
peaks due to finite dot-lead tunnel couplings. I considered both CUE and GOE 
random matrix ensembles of 2D quantum dots. The electron-electron interactions in 
the QD were described by the constant exchange interaction model. S = 0, |, and 
1 spin states of the QD were accounted for, thus, limiting the applicability of my 
results to not so large exchange interaction constants: Jg < 0.5 A. I also assumed 
that T < A < Ec. 

I calculated the ensemble averaged correction to the Coulomb blockade peak 
spacing in the even valley, Eq. (3.41). At Js = my result coincide with that in 
Ref. [43]. I found that the averaged correction decreases monotonically (nonetheless, 
staying positive) as the exchange interaction constant Js is increased, see Figs. 3.5, 
3.6, and 3.7. These dependences are very robust with respect to the choice of RMT 
ensemble or change in the following parameters: charging energy, mean level spacing, 
or temperature. The averaged correction to the odd spacing is of the same magnitude 
and opposite sign. 

I calculated the rms of the correction fiuctuations, see Eqs. (3.57)-(3.64). It is 
of the same order as the average value of the correction, see Figs. 3.5 and 3.6, and 
weakly depend on Jg, see Fig. 3.9. For a small subset of the quantum dot ensemble 
realizations, at the realistic value of Js = 0.3 A, the correction is of the opposite sign, 
see Figs. 3.5 and 3.6. The rms of the correction fiuctuations in the odd valley is the 
same as that in the even one. 

In the experiment by Jeong and co-workers [45] the corrections to the even and 
odd peak spacings due to finite dot-lead tunnel couplings were measured. It was 
found that the even (odd) peak spacing increases (decreases) as the tunnel couplings 
are increased. This is in the qualitative agreement with the theory, see Eq. (3.41). 
The magnitude of the effect was measured at different values of the gas parameter r^ 
(and, hence Jg) as well. Unfortunately, from these measurements one can not draw 
any conclusion about the behavior of the correction to the peak spacing as a function 
of Jg. Hence, more experiments are needed to confirm my theory. To measure the 
ensemble averaged correction to the peak spacing and its fiuctuations as a function 
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of Js one should fabricate a set of 2D lateral quantum dots with different values of 
Ts- Then carry out the measurements on each device under similar conditions. 

In Chapter 5 I studied the effect of mesoscopic fluctuations on the magnitude of 
errors that can occur in exchange operations on quantum dot spin-qubits. I considered 
mid-size double quantum dots, with an odd number of electrons in the range of a few 
tens in each dot. My results indicated two different scenarios for quantum dot qubit 
implementations. 

First, if one is willing to characterize each quantum dot pair separately and have 
them operate one by one, mesoscopic fluctuations will be irrelevant. In this scenario, 
each pair of QDs will require a different pulse shape and duration since QDs are not 
microscopically identical. Multi-electron QDs are tunable enough, easy to couple, 
and much easier to fabricate than one-electron dots; therefore, multi-electron QDs 
are most appropriate for this case. 

Second, if the goal is to achieve genuine scalabihty, one has to operate qubits in 
a similar and uniform way, utilizing a single pulse source. In this case, pulse duration 
time T and switching time r [Eq. (5.9)] should be the same for all QD pairs. Thus, 
four tuning parameters per QD may be necessary to achieve the following goals: (i) 
find isolated, single-occupied energy level (two parameters); (ii) align this level with 
the corresponding level in an adjacent QD (one parameter); (iii) control the inter-dot 
coupling (one parameter). For parameters involved in (i) and (ii), an accuracy of a few 
percent will likely be required. Finally, control over the inter-dot coupling parameter, 
(iii), must allow for the application of smooth pulse shapes in the picosecond range. 
Although these requirements seem quite stringent, recent experiments indicate that 
they could be met [89]. 

In summary, my analysis indicate that mid-size QDs, with ten to a few tens of 
electrons, while not allowing for extremely fast gates, are still good candidates for 
spin-qubits. They offer the advantage of being simpler to fabricate and manipulate, 
but at the same time require accurate, simultaneous control of several parameters. 
Errors related to detuning and sample-to-sample fluctuations can be large, but can be 
kept a secondary concern with respect to dephasing effects provided that a sufficient 
number of independent electrodes or tuning parameters exists. 

In Chapter 7 I studied decoherence of a quantum dot charge qubit due to 
coupling to piezoelectric acoustic phonons in the Born-Markov approximation. I 
calculated the dependence of the Q-factor on lattice temperature, quantum dot size, 
and interdot coupling. The analytical expression for the Q-factor given in Eq. (7.26) 
was found using an expression for the phonon spectral density, Eq. (7.15), which 
takes into account important information concerning the geometry of the double 
quantum dot system. The main conclusion is that, under realistic conditions, phonon 
decoherence is one to two orders of magnitude weaker than expected [81, 104, 105, 
106]. Based on these findings I conclude that phonon decoherence is too weak to 
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explain the damping of the charge oscillations seen in recent experiments [88, 108, 
109]. Hence, other decoherence mechanisms, such as background charge fluctuations 
or noise in the gate voltages, play the dominant role [113]. 

There are two distinct ways to operate a double quantum dot charge qubit: (i) 
by tunnel pulsing or (ii) by bias pulsing. Tunnel pulsing seems advantageous due to 
the smaller number of possible decoherence channels. In addition, the bias pulsing 
scheme, in contrast to tunnel pulsing, introduces significant loss of visibility in the 
charge oscillations. 

Finally, it is worth mentioning that some extra insight would be gained by mea- 
suring the Q-factor as a function of the tunneling amplitude Vm experimentally. Such 
a measurement would allow one to map the spectral density of the boson modes re- 
sponsible for the decoherence. This would provide very valuable information about 
the leading decoherence mechanisms in double quantum dot systems. 



8.3 Possible Directions of Future Research 

In my opinion, the theory of the Coulomb blockade phenomena in the quantum 
dots is in very good shape. I believe there is no experimental result in the last 5 
years which comes as a surprise. Moreover, in most instances there is a quantitative 
agreement between experimental results and the theory. Therefore, one can say that 
the Coulomb blockade phenomenon provide an important tool for analyzing quantum 
mechanical properties of the quantum dot. 

The correction to the spacing between Coulomb blockade peaks due to dot-lead 
tunnel couplings turned out to be very small - both the value averaged over meso- 
scopic fluctuations and its rms. I expect that it is going to be quite challenging to 
discern this effect in the experiment. 

Mesoscopic fluctuations are proved to be an important limiting factor for imple- 
mentation of the quantum dot spin-qubits. However, it should be possible to demon- 
strate basic spin swapping between two quantum dots in the laboratory provided that 
enough independent electrodes (tuning parameters) exist. Besides, the commercial 
pulse generating technique should be signiflcantly improved so that shorter pulses 
with the well-controlled area could be generated. The experiments on the exchange 
operations in the quantum dot spin-qubits are currently on the way. The industrial 
applications require the accuracy of 10^^ and, therefore, are very much questionable 
at this point. 

I find that phonon decoherence is too weak to explain the damping of the charge 
oscillations seen in recent experiments. Hence, more theoretical work needs to be done 
to quantify other decoherence mechanisms. In my opinion, the following mechanisms 
can play the dominant role: (i) background charge fluctuations and/or (ii) noise in the 



105 



gate voltages. On the experimental side, one should measure the decoherence rate as 
a function of the tunneling amplitude. Such a measurement would allow one to map 
the spectral density of the boson modes responsible for the decoherence. In general, 
more work needs to be done to reduce charge decoherence in double quantum dot 
setups. 

In conclusion, I think the field of quantum computing will continue to stimulate 
both the theoretical research and advancements in experimental techniques for years 
to come. 



Appendix A 

Occupation Numbers in the 
Canonical Ensemble 



In the canonical ensemble, where the total number of particles, A^e, is fixed, it is 
difficult to calculate fermionic occupation numbers, Eq. (2.8), directly. This is true 
even in the case of equidistant energy levels separation. 

Fortunately, the bosonization technique allows one to express fermionic field an- 
nihilation and creation operators, ip and ip^ in terms of the bosonic annihilation and 
creation, a^'s and a^'s, and ladder, U and W , operators [34]: 

tlj\x) = ^e-*^^"e-''^'(")t/te-''^("\ (A.l) 

where 

xHx) = ^N + ^J2^me-^^^al; (A.2) 

L is the length of the artificial system; kp is the Fermi wave vector; N = Ne — Ni is 
the number of excess electrons operator; q is the wave vector; and x is the coordinate. 
Since these bosons naturally exist in the grand canonical ensemble, one can fix total 
number of excess electrons, A^, and calculate occupation numbers in the bosonic basis, 
Eq. (2.9). 

Occupation numbers, Eq. (2.9), have the following properties: 

nj = l- rii^j, n^e^^l'^ = n_je~^^/^. (A.3) 

They are similar to np(E) = 1 — nF{—E) and nF{E)e^^^'^ = ni?(— -E)e^^^/^ ones 
of the Fermi-Dirac distribution and refiect the electron-hole symmetry. Combining 
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these two properties one can get the recursion relation: 

Uj+i = 1 - e^^rij. (A.4) 

In the hmit of low temperature, 6E ^ T, for j > we obtain 

n, = e-^' - e-^'^+'^' + O [e-^'^+''>'] ; (A.5) 

in the high temperature limit, T ^ SE: 

1 S sinh [{j - i) f] 



n 



J I'.- 1 



MM + 1 8cosh^[(j-i)|] 



+ O {6') ; (A.6) 



thus, we find first correction to the Fermi-Dirac distribution. The form of the cor- 
rection remains valid even for the slightly non-equidistant energy levels in the dot. If 
this is the case we need to define ratio 6 in the last expression as 5 = SE/T, where 
SE is the mean level spacing. 



Appendix B 

On the Electron-Phonon Coupling 
in GaAs quantum dots 

The material in this appendix is largely based on Ref. [117]. 

B.l Introduction 

There are two types of phonon modes in GaAs/AlGaAs heterostructures that we 
need to worry about: deformation potential and piezoelectric. The former is always 
present in all semiconductor materials, while the latter is due to the lack of inversion 
symmetry in the zinc-blend lattice structure of GaAs. Optical phonons do not usually 
couple strongly to electrons at low temperatures (with some exceptions, like in the 
Peierls stability problem). Acoustic phonons, particularly longitudinal ones, do couple 
more effectively to electrons at low temperatures (transversal phonons do not usually 
couple to electrons in the absence of Umklapp processes). 

B.2 Deformation Potential Contribution to the 
Electron-Phonon Coupling 

The deformation potential only makes sense in the long wavelength, acoustic 
branch limit (see Section 4.12 in Ref. [118] or Section 6.14 in Ref. [119]), where the 
phonons are associated to compression waves of the ion displacement field. Let u(r, t) 
be the displacement field in the continuum limit (so, the crystal is considered as a 
continuum elastic medium). Then, the relative volume change due to the compres- 
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sion wave is given by bV jV = A(r,t) = V ■ u. This volume compression causes a 
local change in the lattice constant, which, in turn, leads to a local variation of the 
conduction band lower edge, Ec'. 

6E,= ^^6V = V^A = AV.u, (B.l) 

where A is called the deformation potential constant (the upper edge of the valence 
band would shift in the opposite direction). Thus, for a plane wave u = uqc**^'^'""'^*-', 
we find that 

5Ec = iAiu-q). (B.2) 

So, it is clear now that only longitudinal deformation potential phonons are important. 
We can put that into a more convenient form for the electron-phonon Hamiltonian: 

i/et-ph = E^^-^(^^)- (B.3) 

i 

We can rewrite the displacement field in terms of normal modes (phonons) and in- 
troduce bosonic creation and annihilation operators: 

u(r) = 5^Qqeqe^''-, (B.4) 

q 
with 

Q^=\hT7^ f^tq + O, (B.5) 



where pion is the ion mass density, l^ is a normalization volume, and Sq is the polar- 
ization unit vector. Also, we can go from first to the second quantized representation 
for the electrons by recalling that, for the Bloch momentum (neglecting Umklapp 
processes) 

{k\e'^-'\k') = Vk'+q. (B.6) 

As a result, 

i^et-ph = E <q (^-q + ^^) 4+q^'^' (B.7) 

k,q 

where the electron-phonon coupling matrix element is given by 



Obviously, this is an oversimplification of the problem since the deformation potential 
constant is actually a tensor in any non-cubic crystal. But since we are only inter- 
ested in the longitudinal modes and have assumed the medium to be approximately 
isotropic, the simplification should be good enough. Notice that |M'^'^^p oc q when 
ujq = sq. Thus the coupling between electrons and deformation potential phonons is 
suppressed at long wavelengths. This is valid, for example, at low temperatures. 



no 

B.3 Piezoelectric Contribution to the 
Electron-Phonon Coupling 

In III-V compounds, there is an electric polarization field associated with the 
strain (displacement) field: 

3 

-Pi = ^ PijkWjk, (B.9) 

j,k=i 

where f3ijk is called the piezoelectric tensor and Wjk is the strain tensor, 



2 \ dxj dxi ^ 

In order to find the associate electric potential energy, we need to solve Poisson's 
equation: 

i,j=l -^ 1=1 i,j,k=l ■' 

Recalling Eq. (B.4) and assuming a normal mode decomposition of the electric po- 
tential, namely, 

^ = Y,^^e^^-\ (B.12) 

q 

we get 

<^q = 4vr^^^3 . (B.13) 

Following the same steps as in the previous section, we can then write 



TT-piezo 
-"el-ph 



J2Kr{bl, + K)cUc^, (B.14) 



k,q 

where the piezoelectric matrix element is given by 

where we have assumed an isotropic dielectric function (e^ = Sije). Now, notice 
that liVfP"'™!^ oc q~^ for tUq = sq. Therefore, the coupling between electrons and 
piezoelectric phonons gains importance at low temperatures. Also notice that M^"'^" 
is real while M'^^^ is imaginary. Therefore, these two electron-phonon amplitudes are 
out of phase by 7r/2 and do not interfere. 



Ill 
B.4 Combining Two Mechanisms 

While the piezoelectric coupling constant is intrinsically tensorial, we can try to 
use an isotropic approximation in order to estimate its magnitude. We introduce the 
substitution 

\ i,j,k=l / q 

where (• • • )q denotes the average over all orientations of q at fixed q = |q|. Thus, for 
the longitudinal acoustic phonons, one can write the total electron-phonon coupling 

as 

This expression should be used with caution since for the averaged piezoelectric con- 
tribution it actually mixes longitudinal and transversal phonon modes. Bruus and 
coworkers [111] go into some details when explaining how the parameter 6 (called 
the piezoelectric coupling) should actually be calculated in terms of phonon velocities 
and the dominant component of the piezoelectric tensor. 

It is convenient to rewrite the electron-phonon coupling Hamiltonian in terms of 
the electron density operator in the momentum representation: 

ffel-ph = J] ^q (bl^ + &q) p(q), (B.18) 

q 

where 

P(q)=E4+q^k. (B.19) 

q 

According to experimental data compiled by Bruus and coworkers [1 1 1] for GaAs 

A ^ 2.2 X 10"^^ J and ^ 5.4 x 10"^° J^m'l 

Thus, one can see that the crossover phonon wavelength, where the piezoelectric 
contribution begins to dominate over the deformation potential, is equal to 

Ac = — = 27r^ ^ 56 nm. (B.20) 



If we recall that s ^ 5 x 10'^ m/s for the longitudinal sound waves in GaAs, we find 
that Ac corresponds to a temperature of Tc ~ 27 K. 

Thus, in the milli-Kelvin range, one can safely neglect the deformation potential 
contribution. 



Appendix C 

Derivation of Eqs. (7.18)-(7.20) 



For t > 0, Eq. (7.2) is time-independent: Hs = Vm^x- Since the K matrix is also 
time-independent [Eq. (7.7)], the matrix A defined by Eq. (7.9) is time-independent 
as well. After some straightforward operator algebra, we find that 

A = - f rfrfi(r)e-^™'""^a^e*"^'""^ (C.l) 



2 JO 



1 f°° 

= - dr B{t) [a^ cos(2t;„r) - ay sin(2t;„r)] . (C.2) 

2 Jo 

One can rewrite Eq. (C.2) as follows 

^ = 2 (^1 + ^^3)^7^ - 2(72 + ^74)0"?/, (C.3) 

where {7i}'s are real coefficients: 

/■oo 

71 + ^73=/ rfr5(r)cos(2t;„,r), (C.4) 

Jo 

72 + ^74= / dTB{T)sm.{2vmT). (C.5) 



The density matrix p{t) is a 2 x 2 Hermitian matrix with unit trace. Hence, it 
has three real independent components and can be written as follows: 

p = - + cr^Repi2 - (Jy\m.pi2 + (yz{pii - -)• (C.6) 

Let us substitute Eqs. (C.6) and (C.3) into the Redfield equation [Eq. (7.8)] and 
use that Hs = Vm^x and K = ^az- A simple algebraic manipulation leads to three 
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differential equations, 








Repi2 = 


= -7iRepi2 + y, 




(C.7) 


Pii = 


= -2vmlmpi2, 




(C.8) 


Im pi2 = 


= (2t;^ + 72)(pii - -) - 


-7ilmpi2. 


(C.9) 



The initial conditions are pii(O) = 1 and pi2(0) = 0. Eq. (C.7) decouples from 
Eqs. (C.8) and (C.9). Its solution is given by Eq. (7.19), where we used the fol- 
lowing identity: 74/71 = — tanh(t>m/T). Eqs. (C.8) and (C.9) form a closed system. 
Their solution is given by Eqs. (7.18) and (7.20). 

The coefficients 71 and 72 [Eqs. (7.22) and (7.23), respectively] are calculated 
using Eqs. (C.4), (C.5), and (7.10). 
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